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Abstract In this second part of the treatment of instantons in quantum mechanics, the focus is on 
specific calculations related to a number of quantum mechanical potentials with degenerate minima. 
We calculate the leading multi-instanton constributions to the partition function, using the formal- 
ism introduced in the first part of the treatise [J. Zinn-Justin and U. D. Jentschura, e-print quant- 
ph/0501136]. The following potentials are considered: (i) asymmetric potentials with degenerate 
minima, (ii) the periodic cosine potential, (iii) anharmonic oscillators with radial symmetry, and (iv) 
a specific potential which bears an analogy with the Fokker-Planck equation. The latter potential has 
the peculiar property that the perturbation series for the ground-state energy vanishes to all orders 
and is thus formally convergent (the ground-state energy, however, is nonzero and positive). For the 
potentials (ii), (iii), and (iv), we calculate the perturbative B-function as well as the instanton A- 
function to fourth order in g. We also consider the double-well potential in detail, and present some 
higher-order analytic as well as numerical calculations to verify explicitly the related conjectures 
up to the order of three instantons. Strategies analogous to those outlined here could result in new 
conjectures for problems where our present understanding is more limited. 
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We continue here the investigations [1] on multi-instantons and exact results. Sections, equations, tables and figures 
are numbered consecutively after those of [1], while the numbering of bibliographic items is independent of [1]. 
For relevant definitions and conventions, the reader is referred to chapters 1 and 2, which are contained in the 
the first part of the treatise [1]. We will often refer to the conjectures summarized in chapter 2.3, and start the 
discussion here with chapter 5. Many of the instanton calculations rely on generaUzations of the ideas and methods 
introduced in chapter 4 (which is again part of the first paper [1] of this series). 
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Chapter 5 



Instantons in General Potentials with 
Degenerate Minima 



5.1 Orientation 

We now consider a general analytic potential possessing two degenerate minima located at the origin and another 
point go > 0: 

V{q) = i(z2+0(g3), (5.1a) 
V{q) = lu>^ {q-qof + o({q-qof) . (5.1b) 
For definiteness we assume u; > 1. 

In such a situation the classical equations of motion have instanton solutions connecting the two minima of the 
potential. However, there is no ground state degeneracy beyond the classical limit. Therefore, the one-instanton 
solution does not contribute anymore to the path integral. Only periodic classical paths are relevant, the leading 
contribution coming now from the two-instanton configuration [2] in the sense of the discussion in chapter 4.4.1. 
Indeed, in the asymmetric case u) > 1, it is useful to redefine the instanton order n in the sense that an instanton 
configuration of order n describes a trajectory in which the quantal particle returns to the original minimum n 
times. 

To calculate the potential between instantons and the normahzation of the path integral, it is convenient to first 
calculate the contribution at /? finite of a trajectory described n times and take the large /? limit of this expression. 
One finds 



y n(l + w) V 2n(l +w) 
with the definitions [see also equations (2.73) and (F.15)] 

f 2oj 



exp I , I exp I -^1^ I , (5.2) 



a = 52^2/(1+-) exp 

Also, 



1 +0; 




y/2V{q) q w(go - q) 



(5.3) 



Aifi)=2 V^Vi<l) dg - 2 ^ ' e-(^/") "/(^+"^ + ■ ■ ■ . (5.4) 
Jo ^ 

Note that n has not the same meaning here as in chapter 4. Since ui is different from 1, the partition function 
may be described, in the path-integral formalism, by trajectories that return to the starting point. Thus, the "one- 
instanton" — ^in the sense of chapter 4.4.1 — configuration does not contribute to the path integral, and n instead 
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CHAPTER 5: INSTANTONS IN GENERAL POTENTIALS WITH DEGENERATE MINIMA 

5.2 The n-Instanton Action 



counts the number of instanton anti-instanton pairs in the language of chapter 4. That is to say, the natural "one- 
instanton" configuration in the asymmetric potential corresponds to the "two-instanton" trajectory in the conven- 
tions of chapter 4. Therefore, n here actually corresponds to 2n in the limit a; ^ 1. 

5.2 The n-Instanton Action 

We now call Oi the successive amounts of time the classical trajectory spends near qo, and tpi near the origin. The 
n-instanton action then takes the form 

n 

A(0i,^j) = na-2^(Cic-"^'+C2e-'^*) , (5.5) 

with Y^"^i {Oi + ifi) = p and [in agreement with (2.60)] 

a = 2 V^V{q) dq. (5.6) 

^0 

By comparing the value of the action at the saddle point, 

n{l+uj) n{l+uj) 
with the expression (5.4), one verifies that one can choose [see also (2.64)] 

Ci = a , C2 = C^/uj , (5.8) 

by adjusting the definitions of 9 and ip. 



5.3 The n-Instanton Contribution to the Partition Function 



The n-instanton contribution to the partition function, in the case of asymmetric wells, has the form 

^ Q-na/g 



(n) 



/3e 



p,—na/g r ( \ 



X exp 



" 1 1 



(5.9) 



The additional term ^ . i (1 — w) 0^ in the integrand comes from the determinant generated by the Gaussian inte- 
gration around the classical path. The normaUzation can be obtained by performing a steepest descent integration 
over the variables Qi and and comparing the result with expression (5.2). The result is 



The factor 1 jn comes from the synnmetry of the action under cyclic permutations of the 0% and y,. 
We now set 



2C 



2C 



A= Ca,\/w, Ml = /"2 = 

T^g 9 ^9 



in such a way that 



-Kg 



(5.10) 



(5.11) 



(5.12) 
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CHAPTER 5: INSTANTONS IN GENERAL POTENTIALS WITH DEGENERATE MINIMA 

5.3 The n-Instanton Contribution to the Partition Function 



As in chapter 4, we introduce the Laplace transform G^"'\E) of {Tre It involves now two integrals, 

I{E - 5, Ml) and 



r+00 

Uy/JI^ / exp{[i(l-a;) + E-i]^-M2e-'^^}d^~J(i;/a;-l/2,M2), (5.13) 
^0 



where I{s, fj,) has been defined in (4.47) and evaluated for — > 00 [equation (4.53)]: 

I{s,ii)^li'+'/^T{-s). (5.14) 

Then, G^") {E) can be written as 

G(")(£) = iA" 1,^2)]" • (5.15) 

From the sum over n, one then infers A(£'), the sum of the leading multi-instanton contributions to the Fredholm 
determinant ^(i^). One finds 

Aw= ,„ ' ^ .(-?^Y(-?^Y"t;!i, (5.16, 



r{^-E/cj)T{^-E) \ g J \ cog J 27r 

where now we have added the harmonic oscillator contributions corresponding to the two wells. A generalization 
of this equation, including effects of higher order in g, leads to (2.64). 

This equation A{E) = has two series of energy eigenvalues, close for 5 — > to the poles 

En = N+l + 0{X), (5.17a) 
En = (TV + i) a; + C(A) , (5.17b) 

of the two F-functions. The same expression contains the instanton contributions to the two different sets of eigen- 
values. One verifies that multi-instanton contributions are singular for a; = 1. However, if one directiy sets w = 1 
in the expression (5.16), one obtains 

1 / 2G\^^ p-a/g 

an expression consistent with the product A+(£) A_ (E) of the results (4.60) found for the double-well potential, 
and hkewise consistent with (2.21). 
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Chapter 6 

Instantons in the Periodic Cosine Potential 



6.1 Orientation 

Analytic periodic potentials lead to additional complications, which the example of the cosine potential illustrates. 
We consider the Hamiltonian 

9 f d\\ 1 



where the peculiar normalization has been chosen to avoid the proliferation of big integer factors. 
The unitary operator T that translates wave functions by one period T = 7r/2 of the potential, 

Ti;{q) = i;{q + T), (6.2) 

commutes with the Hamiltonian. Because the adjoint of T fulfills T^tp{q) = ijj{q — T), the eigenvalues of T have 
unit modulus. Indeed, T can be diagonalized simultaneously with the Hamiltonian, 

TV'(g)=e'^V(9)- (6.3) 

This corresponds to decomposing the initial Hilbert space Ti into a sum of spaces H^. The coordinate q can then 
be restricted to one period or more conveniently considered as an angular variable parameterizing a circle. 

In Tiif the Hamiltonian has a discrete spectrum with eigenvalues Eff{g, (p) and eigenfunctions V'Af,<^(</)- 

TtjjN,^ = e"^tjjN,^ , Hil}N,ip = EN{g, <f) V'iv.y • (6.4) 

The eigenvalues E^ig, are periodic functions of the angle (fi and Ei^{g, ip) = N + ^id)- 
To each state of the harmonic oscillator is thus associated in H, for g small, a band when <^ varies. 

6.2 The Partition Function in the (/^-Sector 

We now introduce the partition function in H^: 

2(/3,5,<p)=^e-'^^"(s,<P). (6.5) 
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CHAPTER 6: INSTANTONS IN THE PERIODIC COSINE POTENTIAL 

6.3 Perturbation Theory and Instantons 



In the initial Hilbert space H we also define the quantity 

Zi{(3,g) = Tr(T'e-^^) 

dq {q\T'e-^" \q) = I Aq {q + lT\ e-^" \q) 
Jo 



- hi 



Inverting this last relation, one finds 



+00 



(6.6) 



(6.7) 



This is a generahzed partition function with twisted boundary conditions, which depends on the rotation angle ^. 
The path integral representation of Zi (/3, g) can be written as 



J Q\ 



with 



g(/3/2)=g(-/3/2)+;7r/2 
(3/2 



[dg(i)]exp[-5(g)/5] 



At 



13/2 



^ + 3^ (1- COS 4?) 



(6.8) 



(6.9) 



The integration variable q parameterizes a circle and the boundary condition is periodic with the constraint that 
q{t) belongs to the topological sector of trajectories turning I times around the circle. 

Note that since 

21 2 i /■+'^/^ 

U^=-^ {q{P/2) - q{-P/2)) = - if dt q{t) , (6.10) 

7^ 7^ J-0/2 

the factor exp[i I can be incorporated into the path integral by adding a topological term, the integral of a local 
density, to the action: 

9 9 J-B/2 



(6.11) 



The sum (6.7) can thus be written as 

9,f)= j [dq{t)] exp 



5(g) , 21 



+13/2 



P/2 



dtq(t) 



(6.12) 



with now unrestricted periodic boundary conditions on the circle: q{(3/2) = q{—(3/2) (mod 7r/2). 



6.3 Perturbation Theory and Instantons 

In an expansion for /3 ^ oo, due to the boundary conditions in the expression (6.9), Zo{0, g) is dominated by the 
perturbatlve expansion and Zi{(3, g) for Z 7^ by Instantons: In a band all eigenvalues are degenerate to all orders 
in perturbation theory. Only for i = ±1 do instantons correspond to solutions of the classical equation of motion: 

qc{t) = arctane='=* S{qc) = | . (6.13) 
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CHAPTER 6: INSTANTONS IN THE PERIODIC COSINE POTENTIAL 



6.4 Multi-Instantons 



For |Z| > l,Zi (/3, g) is dominated by multi-instantons. 

The energy eigenvalue {g, can be expanded in a Fourier series: 

+00 

Ejv(5,¥')= 1] i^JV,;(5)e'"^, En,i = En,-i . (6.14) 
/=— 00 

For g small, E^^i {g) is dominated by Z-instanton contributions. In particular, for the ground state energy Eq (<?, f) 
in the if sector, the I = 1 term behaves like 

Eo,i=i{g) ~ ^e-V2s. (6.15) 



6.4 Multi-Instantons 

There is one important difference between the double well and the cosine potentials. In the case of the double 
well potential, each configuration is a succession of alternatively instantons and anti-instantons. Here, by contrast, 
the paths consist in an arbitrary succession of turns in the positive and the negative direction, that is an arbitrary 
succession of instantons and anti-instantons. Therefore, we assign a sign e = +1 to an instanton and a sign 
£ = — 1 to an anti-instanton. A straightforward calculation, similar to the calculation presented above (for details 
see Appendix C.2) yields the following interaction term between two consecutive instantons of types £1 and £2 
separated by a distance 612: 

^e-^-. (6.16) 
9 

The interaction between instantons is repulsive, while it is attractive between instantons and anti-instantons. 

We redefine the parameters A and jj. previously introduced in equation (4.42) in the context of the double-well 

problem, 

A=^e-i/29, ^=-. (6.17) 



/27r g 
The one-instanton contribution at leading order can then be written as 



/ng 

With this notation, the n-instanton contribution reads 



^c-i/25 = AVm- (6.18) 



Z^^\f3,g,ip)=Pe-^/^— [ (5 ( V^i-/?) JnW (6.19) 



with 



Jn{G)= Yl exp ( ^--£i£i+ie-^- +i£iV3 J . (6.20) 

The additional tenaisup comes from the topological term in the expression (6.12). We have identified £„+i and 
£1. 

Since the interaction between instantons contains both attractive and repulsive contributions, we sum series with g 
imaginary and then perform the analytic continuation of both the Borel sums and the instanton contributions. 

Following the same steps as in the case of the double-well potential, we introduce the contribution G'-"^ {E, g, tp) to 
the resolvent, Laplace transform of 2^^"^ (/?, g, (p). The integral over the 0i again involves only the function (4.47), 
which we replace by its asymptotic form (4.53). We obtain 

n dE 



G^-\E,g,^) = {[T{l-E)T„ 



X ^ exp 

{£i=±l} 



Y'lSiip + {E - \) \n{ei£i+i) 



(6.21) 
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CHAPTER 6: INSTANTONS IN THE PERIODIC COSINE POTENTIAL 



6.4 Multi-Instantons 



where the choice in the determination of ln(eiei_|_i) depends on the initial phase of g. We choose 

ln(eiei+i) = -iiTr (1 - eie^+l) . 
The expression (6.21) can then be rewritten as 



(6.22) 



X ^ exp 

{£i=±i} 



^\£i ip - \ \TT {E - \) (I - £i£.,+ i) 



(6.23) 



The summation over the set {si} corresponds to calculating the partition function of a one-dimensional Ising model 
with the transfer matrix 



M 



-i7r(B-l/2) 



-i7r(E-l/2) 



The sum then is simply Tr M". The expression (6.23) then becomes 

G'(")(E,5,^) ~ — ^ {[r(l -E)/x^]"TrM"} . 



(6.24) 



(6.25) 



The sum 0{E,g,(p) of all leading order multi-instanton contributions can now be calculated. One finds (using 
Indet = Tr In) 

g{E,g,^) = -^lndet[l-Xr{^-E)t,''M] . 
The contribution to the Fredhohn determinant (fi) in the v'-sector thus is 
A{E, if) = det[l - XT{1 - E) i^E - M] 

= l-2Acos^r(i-E)M^ + A2 [Ti^-E)]\^^ + . 
Again, after addition of the partition function of the harmonic oscillator, the expression becomes 



(6.26) 



(6.27) 



A{E,ip) = 



i 2 cos (fi / 2\ i/2g 



TU-E) \g 



+ 



-1/2S / _ 



-l/2s 



9 J ni+E)' 



(6.28) 



Remarkably enough, the equation A{E, y>) = can also be written in a form that is symmetric in the exchange 

g,E^ -g,-E: 



.1/23 



+ - 



-1/2S 



2 cosy' 

gj T{^-E) ' V gJ r{^+E)~ ^ ' 



(6.29) 



a property that depends expUcitly on the normaUzation of the one-instanton contribution. This symmetry, however, 

is slightly fictitious because the equation is actually quadratic in r(i — _E) and only one root is relevant for g > 0. 
The spectral condition (6.29) is fully compatible with (2.80) at leading order in g. Indeed, the generalization to 
higher orders in g is 



2X-Bpc(E,s) gApc(B,s)/2 / 2\Bi-(S,s) g-Apc(B,9)/2 2 COS 

r[i + Bpe(£;,5)] " "TIT 



r[i-Bpe(£,ff)]^ V 5 

We supply here a few more terms in comparison to (2.81), 

1 



(6.30) 



B^,{E,g) = E^\E 



25 . 35 



2 

721 



■E' 



17 



4 32 

,4 , //i/-„5 



(6.31) 
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CHAPTER 6: INSTANTONS IN THE PERIODIC COSINE POTENTIAL 



6.4 Multi-Instantons 



The first few terms of the expansion of the instanton ^pc-function are: 

A^{E,g,j) = ^+(3E^ + ^^ g 



„ 215 

^ H 

64 



+ (11^3^1^) ^,^^_199^,^341^, 
+ (^^E' + ^2&E' + ^E^g^ + 0{g') . (6.32) 
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Chapter 7 

Instantons in Radially Symmetric 
Oscillators 

7.1 Orientation 

It is interesting to consider a last example, the analytic continuation of the energy eigenvalues of the 0{u) sym- 
metric anharmonic oscillator corresponding to the Hanoiltonian 

if = -iV2 + iq2+5(q2)^ (7.1) 
from g > to g < 0. The radial one-dimensional Hamiltonian is given in (2.84), 

rr, s lu-ld ll(l + iy-2) 1 2 4 

For g > 0, the potential is bound from below, and there are no degenerate minima. Therefore, the quantization 
condition for g > 0, which is given in (2. 101), can be considered as rather trivial: only the perturbative i?-function 
enters, and since there are no degenerate minima, there are no nontrivial saddle points of the Euclidean action, and 
thus no instantons to consider. 

For g < 0, there are two possibilities: one may either endow the Hamiltonian with a self-adjoint extension, as 
discussed in chapter 2.2.4 and below in chapter 7.2, or one may consider resonances, as in chapter 2.2.5 and below 
in 7.3. 



7.2 Self-Adjoint Extension 

7.2.1 Double- Well, Symmetry Breaking and 0(z^) -Potentials 

One particular aspect facilitates considerably the analysis of the self-adjoint extension of (7.2): There is a general 
connection [3-6] between the double-well potential with an additional symmetry breaking term ("broken-double- 
well") on the one side and anharmonic oscillators with 0(i^)-symmetry on the other side (acting in IR''). Typically, 
the sign convention in this investigation is chosen such that the double-well potential is bounded from below, and 
the anharmonic (^(i/j-symmetric oscillator is formulated for g < 0. We here briefly explain the relevant method 
and recall the radial Hamiltonian (2.84) with negative couphng 



For iUustrational purposes, we have explicitly replaced g ^ —g in comparison to (2.84). The potential, in this 
convention, is not bounded from below for a (redefined) positive g. 
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CHAPTER 7: INSTANTONS IN RADIALLY SYMMETRIC OSCILLATORS 

7.2 Self-Adjoint Extension 



We now try to find a self-adjoint extension of this Hamiltonian; although the coupling is negative, a self-adjoint 
extension implies that the spectram will be real. The Hamiltonian (7.3) leads to a radial Schrodinger equation 



-/'(r) - 1 ^'[r) + ^-^ii^ ^ i,[r) + (r^ ~2gr^) ^j{r) = 2E ijj{r) . 



(7.4) 



In order to show the connection between the "broken-double-well" potential and the anharmonic oscillator, we use 
here a derivation based on differential equations, but the same results can be obtained by path integral methods (for 
an overview of related issues see [4]). We consider first the more general Hamiltonian (3.78) 



d l{l + v-2) 
r dr 



+ - F(r) . 
9 



(7.5) 



We have used here the scaling x ^ xj y^. In the case of (2.84), we would have V{r) = + r^. Because V{r) 
is an even function, we now set 

^2 



V{r) = W^(r2) = — + 0{r^) . 



We first eUminate the 1/r^ term by the change of variables 



and obtain 



-| [X"(r) 
where we have set [see also (2.87)] 



2j + i ,( ; 



-^W{r^)x{r)=Ex{r), 



We then take = a; as a new variable, 
and find 



-g [2 X £,"{x) + 2 (i + 1) (,'{x)] + - W{x) ^{x) = E ^{x) . 



We now write ^ as a Laplace transform 



^(x) = J dp(f"''B'n{p). 

With proper boundary conditions, the equation becomes 



W -g— + 



g d 



dp J 2 dp 



V{P) + 2 (V - 1) ri'{p) - 2 (i - l)p7){p) = Erjip) . 



(7.6) 

(7.7) 

(7.8) 

(7.9) 
(7.10) 
(7.11) 

(7.12) 
(7.13) 



In the special example W{x) = x/2 — x^ given in (7.4), to which we restrict the discussion in the sequel, the 
equation is again a second-order differential equation. After a last transformation to eliminate the r/' term. 



7?(p) = a(p)(^(p), with 



1 - V 
45 



one finally obtains 



2^"(^)+2^ 



V 



E 



(fiip) - j P (fiip) = T7Vip) 



(7.14) 



(7.15) 



This is equivalent to the Schrodinger equation for a particle moving in a double-well potential with minima and 
p ± 1/2, and an additional symmetry-breaking term j p. The equation has been derived by appropriate substitutions 
in (7.4) and naturally provides a self-adjoint extension of the Hamiltonian (7.3). This concludes the derivation of 
the correspondence between the anharmonic oscillator at negative couphng and the double-well potential with a 
linear breaking term, which is expressed by the equations (2.90a) and (2.90b). 
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CHAPTER 7: INSTANTONS IN RADIALLY SYMMETRIC OSCILLATORS 

7.2 Self-Adjoint Extension 



For j = 0, the symmetry-breaking term vanishes. In this case, an comparison of (7.15) with (2.10) leads one auto- 
matically to the A and B functions of the transformed problem (7. 15) . These have to be the same as those for the 
anharmonic oscillator obtained in (2.89a) and (2.89b) as we set j = 0. We thus recover the correspondence 
between the double-well potential on the one side and the anharmonic oscillator with radial symmetry, considered 
in the case of negative coupling and endowed with a self-adjoint extension, on the other side, 

B4E,g,j = 0) = 2Bd^{E/2,-g), A,{E,g,j = 0) = Ad^{E/2,-g). (7.16) 

This has already been mentioned in chapter 2.2.4 [see equation (2.91)]. This correspondence is valid on the per- 
turbative level (B function) as well as on the level of instanton effects (A function). For j ^ 0, the reflection 
symmetry p —p is broken at the same order as the first quantum correction, so that the degeneracy is lifted. 



7.2.2 Degenerate Minima and Symmetry Breaking 



We now consider a slight generalization of the situation encountered in the Fokker-Planck equation (2.2.4), where 
the potential has a symmetric structure broken at relative order g, and may allow for a vanishing perturbative 
expansion. Indeed, in several situations, for instance in the case j — ±1 in (7.15), the potential may have the 
typical form [we assume the general structure (2. 10a) for the Hamiltonian] 



Vtot{q) = V{q)+gVp^rt{q): 



(7.17) 



where the form y is a synnmetric potential [V{q) = V{qo — q)] with degenerate minima and V^ert(9) breaks the 
symmetry at order g. When the breaking term is treated at leading order, it is only the difference between the values 
at the minima of the symmetric potential which is important. Therefore, as a simplifying feature but without loss 
of generality at leading order, we consider the example of a Unear synmietry breaking potential 



Vtot{q) = V{q)+gri- 



qo 



(7.18) 



where we assume V{0) = 0,r] > 0. 

The analysis is, in many respects, analogous to the calculations of chapter 5.3 (see also reference [2]). Eigenvalues 
are only degenerate at leading order, and the degeneracy is Ufted at order g, the ground state corresponding to the 
well at g' = 0. At leading order, the instanton contributions can be derived from chapter 5 with w = 1. Again, we 
call ipi the times spent near q = and 9i near q = q^. Then, the n-instanton action reads 



A{0i,ifj) =na-2C 

with Er=i i^i + <Pi) = and [see also (2.58) and (2.60)] 

rqo 

a = 2 ^/2V{q) dq, 
Jo 



(7.19) 



C = qo exp 
The n-instanton contribution then has the form 



dq 







qo 



(7.20a) 
(7.20b) 



L J (n) 



^^-0/2 



X exp 



n 

-rj^ei--A{9,^) 



(7.21) 



The additional term —rj 9i in the integrand comes from the value of the potential at g = go- 
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The only required integral [equation (4.47)] has already been evaluated. The sum of leading order multi-instanton 
contributions reads 

1 / 2C\^^"'' e-^/s 

^(^)° r(i-^)r(i.,-.) n-T) 

Note that when rj is a positive integer, the one-instanton contribution to the states Ej^ = N + ^is real, indicating 
that the behaviour of the perturbative expansion at large order k is at least smaller by a factor 1/fc than naively 
expected. For < 77, the instanton contribution is of order e~"/^, while for > 77 it is of order e~'^^^^. 



7.2.3 The Fokker-Planck Hamiltonian 



A simple example illustrating the remarks of chapter 7.2.2 is provided by the Fokker-Planck Hamiltonian (see 
chapter 2.2.4): the stationary solution, as shown below, is not normalizable, and instanton effects determine the 
energy of the ground state [10]. 

We recall the Riccati equation (3.14), 

gS'{q)-S\q) + 2V{q)-2gE = 0. (7.23) 
This equation formally allows for a solution with = if the potential V{q) has the following structure, 

V{q)^^[UHq)-gU'iq)] . (7.24) 

Indeed, this structure automatically leads to a class of potentials for which the perturbative expansion of at least 
one eigenvalue vanishes identically to all orders in the coupling. We assume U (q) to be a polynomial that has two 
zeros at g = and q = qo > 0, such that 



U{q)=q + 0{q^), U{qa - q) = U{q), 

and assume that it is an exact solution of the Riccati equation (3.14) with E = Q. 
We now speciaUze the treatment to the Fokker-Planck potential 

2 



(7.25) 



1 



Vfp(p) = 1: [p 



+ 9P, 



(7.26) 



which is obtained in a natural way by setting j = — 1 in (7. 15). In a perturbative expansion around the well p = - 
one finds a ground state with = to all orders. The equation may alternatively be written as (p = g — 1 /2), 



(7.27) 



At leading order in g, V{q) is symmetric with degenerate minima. The breaking term lifts the degeneracy at order 
g and implies that the perturbative ground state, corresponding to the well at g = 0, has E = while the lowest 
state in the other well has E = 1. 



In fact the issue is more compUcated because the wave function 



ip{q) = exp 



1 

--J dq'U{q') 





1 








'1 






= exp 




a 


4)] 


oc exp 










.9 








.9 







(7.28) 



is not normalizable, and thus is not an eigenfunction. An analogy with the Fokker-Planck equation suggests that 
the case E = Ohe identified with an equilibrium probability distribution. Therefore, the non-normalizable wave 
function (7.28) may naturally be identified with a "pseudo-equilibrium" distribution. 

The true ground state has E > and at leading order for 5 ^ is dominated by the one-instanton contribution. 
Shifting E in the expression (7.22) as E ^ E + 1/2 and setting r] = 1, one finds 



A{E) 



T{-E)T{1 - E) 



_^ 2 



2E 



-1/39 



27r 



(7.29) 
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the generalization being (2.96) 



/_2y--(-^^)exp(-A.p(i.,,))^^_ (7.30) 



T{-BME,g))r{l-BFp{E,g)) ' \ gj 2tt 
At leading order, for the ground state, 

e-i/3s 

^0(5) (7.31) 

An imaginary part appears only at two-instanton order and governs the large-order behaviour of the non-Borel 
sununable one-instanton expansion, where for the classification of the instanton order we follow the convention as 
outlined in chapter 5.1. The arguments presented here can be easily generalized to the non- symmetric situation. 

The perturbation series, in the Fokker-Planck potential, vanishes to all orders for the ground state only. For excited 
states, the leading-order as well as corrections of relative order g,g-^,g^, . . . are nonvanishing. The following 
statements illustrate this phenomenon. First, we supplement the perturbative Bpp-function given in (2.97) by an 
expression valid up to the order g^, 

BFp{E,g) = E + 3E^g+(35E^ + ^E^ 5^ 

+ (1^ + 105 E=) »3 + + 1521^ E= + if! ^) (7.32) 

Inverting the perturbative quantization condition (2.97) 

BFPiE,g)=N (7.33) 
for general N, we obtain the following perturbative expansion up to 0{g^), 



EN{g)^N-3N^g- (^17 N^ + ^N^ g^ 



375 . 165 , / 10689 9475 ,^3 1105 ^ A 4 

— N^ + —N^\ g^- f ^_iV^ + ^iv3 + -_iVj g\ (7.34) 

For the ground state (N = 0), all the terms vanish, whereas for excited states with N — 1,2,..., the perturbation 
series is manifestly nonvanishing. For completeness, we also supplement the instanton App-function (2.99b) for 
the Fokker-Planck potential by a few more terms, 

AppiE,g) = j-+(^17E^ + l^ g+(227E'+^-^E^ g^ 
f 47431 . 11485 1105\ 3 

+ + + + (7.35) 

The functions Bpp and App determine the perturbative expansion, and the perturbative expansion about the in- 
stantons, in higher order (see also chapter 8 for an appUcation of related ideas to the double- well problem). 

(Remark.) When investigating a shght generalization of the potential (7.26), namely [see also (2.90b)] 

yFp(p) = l {p^-\) +9JP, (7.36) 

it may be shown that for negative integer j, the perturbation series of at least one energy level terminates. Assuming 
that the solution of equation (7.11) is regular at a; = 0, we expand it, as well as the potential, in a Taylor series, 
after changing variables x 1— > gx, 

l^(x) = ^lK<?"x", ^(x) =^en<?"a;". (7.37) 

n=l n=0 



18 



CHAPTER 7: INSTANTONS IN RADIALLY SYMMETRIC OSCILLATORS 

7.2 Self-Adjoint Extension 



We then obtain the recursion relations 

n 

-2g{n+ 1) (j + 1 + n) U+i+Y, 5""^"' Wn-p ^p = EU. (7.38) 

If j is a negative integer: j = —N — 1 (with > 0), the set of linear equations with < n < is closed. Idem 
est, the coefficient j + l + n vanishes in this case for n = N, and = for all n > A''. Alternatively, one observes 
that the determinant of the [A^ + 1] x [A^ + 1] system of linear equations defined by equation (7.38) has to vanish, 
a condition that determines -E as a solution of an algebraic equation of degree A^ — 1. This property reflects the 
equivalence of 0{i') models with v even and negative with fermion Hamiltonians. For N = and thus j = —1, 
one recovers the known result 

E = 0, (7.39) 

which holds in a strict sense only on the level of perturbation theory [see also equation (7.31)]. For j = —2, one 
finds E = ±1. The instanton corrections follow directly from equations (7.22) and (7.29). As an illustration, at 
next order, and thus, for j = — 3, one finds 

E^ -AE-B,g W"{Q) = , (7.40) 

an equation that has three solutions for small g. Also, for g = Q, the solutions are = 1 — |j | + 2fc with k integer, 
< < b l — 1- We conclude this chapter by noting that the Hamiltonian (7.15) impHes a spectral symmetry with 
respect to the sign change j —j. However, here this symmetry is broken in the discussion following (7.38), 
because the potential (7.26) has a unique global minimum at p = —1/2, not p = +1/2 [analogous statements 
hold for the generalization (7.36)]. For this reason, the system of equations is closed for negative integer j (but the 
closure does not hold for positive j). 

For a more detailed treatment of the intriguing aspects related to the Fokker-Flanck potential, the reader is referred 
to [11]. 

7.2.4 0{i') -Symmetric Quartic Potentials 

We have related in chapter 7.2.1 perturbative expansions of 0{u) synnmetric quartic potentials and double-well 
potentials with Hnear symmetry breaking. The additional contribution to the potential — jp has the effect of adding 
a contribution ±j/2 to the action depending whether the instanton is close to 1/2 or —1/2, respectively [12]. 
We now call 6i the successive amounts of time the classical trajectory spends near —1/2, and tpi near 1/2. The 
n-instanton action then takes the form 



A{9i,^,) = ^~2^{e-'^+e-'^') (7.41) 
i=i 

with Er=i i^i + Vi) = /3- We set 

e-1/35 2 
A=^— , M = — • (7.42) 
27r g 

The n-instanton contribution then takes the form 



X exp 



Y^-j{ei-^i)--A{e,ip) 

.i=l ^ 



(7.43) 



The sum then involves the integrals (taking into account E t-^ E/2) 

^ e^v{[±\j + \E-\)\e-HG-')de=X{\{±j + E-l),,i). (7.44) 

Jo 
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Using the asymptotic form (4.53), 

/(s,/x)~M^+'/'r(-s), (7.45) 

we obtain the spectral equation 

1 / 2\^e-V3s 

AJE) = —r, ^-n r + — = . (7.46) 

The generalization to higher orders in g is given by (2.92) and reads 

1 

r [i(l + j - B,{E, -g,j))\ r [i(l - J - B,{E, -g,j))\ 

D (IT ^ A\ 



^\~~gj 2^ 



We 

(B.12)] 



supplement here a few terms in comparison to (2.89). The perturbative i?^ -function reads [see also equation 
12)] 



B^{E,g,j) = E+(^-^E^ + ^-^-^^ g 

f35 25 ^ 15 .2 ^\ 2 

+ — H E r E] g^ 

+ I E'' - — E^ + — f E^ -— + — f - —] 5^ 
' 16 8 8^^ 16 •^ 8-^ 16'^ 



1155 ^4 735 ^2 315 35 105 .n 175' 

E^ E^ H f E^ r H r 

^ 16 8 8 le-' 8 16 ^ 

/ 45054 . 45045 31185 „ 15015 , „s 

+ E^ H E-^ H E 7^ E-^ 

V 64 32 64 32 

^?^fE-l^fE),' + 0(,^y (7.48, 
The first few terms of the expansion of the instanton Ajy-function are: 

(227 77 .2 ^ 187 \ ^ 

/ 47431 ^4 _ + 3717 2^2 _ 341 , 

\ 192 96 32 64 

1281 .2 28829\ 3 
32 ~ 576 ) ^ 
/ 317629 264725 . 842909 19215 . ^ 

(7.49, 

In (2.89a) and (2.89b), the functions and Aj, are given up to terms of the order g"^, respectively. 



7.3 Resonances 

7.3.1 The e)(2)-Aiiharmonic Oscillator 

We first discuss resonances in the example case u = 2 (see also chapter 2.2.5), the generalization then being 
simple. Note that for z/ = 2, we have j = 1 + v/2 —1 = 1 [see equations (2.87) and (7.9)]. As in chapter 7.2.2, we 
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extensively use here the ideas outUned in chapter 5.3. For g > 0, a& discussed in chapter 7.1, there are no saddle 
points of the Euclidean action beyond the trivial one, and thus no instantons to consider. 

For < 0, the instanton solution exists and has the form 

q(t)=u/(i), (7.50) 
in which u is a fixed unit vector. The leading one-instanton contribution to the ground state energy is 

Im£;(i)(5) = -e^/^s (l + 0(£f)) forg^O-. (7.51) 

a 

It is easy to calculate the instanton interaction, and thus n-instanton action 

A{9i) = —n- 4 J2e~^' cos ipi, (7.52) 

i 

in which 9i is the distance between two successive instantons and cpi the angle between them: 

cos (pi = Ui- Ui+1 . (7.53) 

It is convenient to consider the quantity 

Z{0, a)=Tr [R{a)e-''"] = J [dq(t)] exp [-S{q)/g] , (7.54) 

where R{a) is a rotation matrix which rotates vectors by an angle a. In the path integral, it leads to the boundary 
condition that q(t) at initial and final times differ by an angle a: q(— |/?) • q(|/?) = cos a. 

We recall that in two dimensions with a radial coordinate p and an angle cf), we have = p~^{d/dp) (p d/dp) + 

{d'^/dcj)'^). Eigenfunctions of the angular part are of the form e"''"^ with I = —oo, . . . ,oo. The r.h.s. of 
equation (7.54) can be rewritten as 

Tr [R{a) e^/S^] = ^ e"' ' ""'^ ^'-^ , (7.55) 

l,N 

where I is the angular momentum. The boundary condition in the path integral (7.54) imphes for the multi-instanton 
configuration the constraint 

n 

^ = a , (7.56) 

i=l 

which can be implemented through the identity 



(n \ +00 " / n 

E^'-" =^ E E 
i—l / 1— — 00 L 1 



ifi -lal 



(7.57) 

As before, we introduce the Laplace transform of the twisted partition function Z{(3, a) defined in (7.54), 

roo 1 

G{E, a)= d/3 e^^ Tr [R{a) e"^^] = E e"" " ^ ■ (7-58) 



l,N 

The n-instanton contribution to expression (7.55) then takes the form 

+ 00 



G(^){E,a)r~.-^ ^ e-''«[/(E-l,Z,M)r, (7.59) 

/— — oo 



where we have again redefined A and fj, [cf . equations (4.42) and (6. 17)], 



A = -ie^/^», /u = --, (7.60) 

9 

I{s,l,ii) = 7^-1 dip dO exp [s 9 + iltp — fie~^ COS If) . (7.61) 
271' Jo Jo 
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We introduce the generating function of n-instanton contributions at fixed angular momentum I: 

SiiE,9)=-^ ln[l + XI{E-l,l,n)] . (7.62) 

To evaluate I{s, I, fi), we first integrate over 6, where we employ the approximations — > oo and thus g — > 0_, as 
in going from (4.51) to (4.52). We find 

7(s,/,M) = M^+'r(-s) _^'"^e''^ {cos^fiY. (7.63) 

The last integration yields 

Jo 2^ ^' r(i(s-j) + i)r(i(s + i) + i) 

and thus 

:i(5-0 + i)r(i 



'^^'''^^ I2; sin(7rs) r(i(s-o + i)r(i(s + o + i) 

r(i(s + o + i)' 



^2^ rri(s + n + ir ^ ^ ^ 



the first expression being explicitly even in /. 

The sum of leading order contributions to the Fredholm determinant at fixed angular momentum / thus is 

The eigenvalues Ei^n are solutions of the equation A; (E) = that satisfy 

Ei,N = l + '2N+l + 0{g), N>0. (7.67) 

7.3.2 The C(i/)-Symmetric Hamiltonian 

One can extend this result to the general O(i') case since, at fixed angular momentum /, the Hamiltonian depends 
only on the combination j = I + 1^/2 — 1, and thus the generalization of the discussion in the previous chapter 
simply involves the replacement I j. Hence, making in equation (7.66) the corresponding substitution, one 
obtains 

This equation is consistent with (2.102) at leading order in g. The imaginary parts of the energy eigenvalues 

Ej,N =j + 2N+l + 0{g) , iV > , (7.69) 

for 5 — > 0_ follow: 

'^"^^'^^o^-mrUTYTW) [-g) ^vs. (i + (,)) , (7.70) 

in full agreement with the ground-state result (7.51). Using the Cauchy formula, one can then derive from this 
expression large-order estimates for perturbative expansions [13]. At next order in A, one obtains the two-instanton 
contribution, which is related by the same dispersion relation to the large-order behaviour of the perturbative 
expansion around one instanton. 

Finally, note that checks about these expressions are provided by the perturbative relation between the 0{i^) an- 
harmonic oscillator with negative coupling and the double- well potential with Unear symmetry breaking derived in 
chapter 7.2.2. 
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Chapter 8 

Instanton Calculations in the Double-Well 
Problem 



8.1 Orientation 

The purpose of the current chapter consists in the illustration of the general discussion of instantons by concrete 
calculational example. We choose the quantum-mechanical double-well oscillator with a potential of the form 
^{q) = 9^ (1 — as in equation (2.10). The resurgent expansion for the energy of a level with principal 

quantum number N and parity e is then given by equation (2.13). We remember the nonperturbative factor (2.27) 



1 

which characterizes the instanton contributions and the logarithm (2.28) 



(8.1) 



which generates, for g positive, imaginary contributions that cancel among the resummed perturbative expansion 
and the instanton contributions. 

The calculations, and even the results which are discussed in the current chapter, have a somewhat involved struc- 
ture and are rather lengthy. The formulas should be taken cum grano salis, exemplifying the application of general 
concepts discussed in previous chapters to a case of special interest. The lengthy formulas and numerical results 
(see also Table F.l below) are not displayed for their own sake; they receive a meaning and an interpretation in the 
context of the general conjectures introduced and motivated in chapters 2, 3 and 4. 



8.2 Corrections to Asymptotics 

According to (2.14), the perturbation series for the level N is independent of the parity e and can be written as 

oo 

<(5)=E<k5^- (8-3) 

We have determined the first 300 terms of the perturbative expansion of the ground state(s) with A'' = (e = ±) 
in closed analytic form (a complete list is available for internet download at [7]). The coefficients are expressed in 
terms of rational numbers. An example is given in equation (2.7). The first terms read [see also equation (2.6)]: 

r.(o)/ N 1 9 2 89 3 5013 4 88251 . 3662169 g 7, ,^ 

Eri9) = ^ - 9 - ^9^ - Y 9 - 9^ - 9' ^ + ' ^ ^ ^ 
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The magnitude of the terms grows factorially, and the series is nonaltemating: all coefficients except the first are 
negative. We had stressed in chapter 2.1.1 that the imaginary part incurred by analytic continuation from g negative 
is compensated by the explicit imaginary part of the two-instanton contribution to the energy eigenvalues. For the 
ground state(s) with N = 0, the formal expansion of the two-instanton effect as given by equation (2.15) reads 



-l/3s 



ng 



In 



K 



eo,2o/f 9 



K=0 



The imaginary part (generated by the logarithm for positive g) reads 



-l/3s 



[eo,210 + 60,211 g + 60,212 5^ + 0{g^ 



(8.5) 



(8.6) 



The leading factorial growth of the perturbative coefficients Ej^]^ for the ground state is known [2] to be of the 

form ^ — S^"*"^ Kl /tt. By analytic continuation of the perturbation series for positive g, an imaginary part of 
the order of exp(— 1 / 3g) is obtained. The large-order behaviour of the perturbative coefficients and the imaginary 
part of the two-instanton energy shift are connected by the (dispersion-type) relations (2.42) — (2.44). Indeed, as 
discussed in chapter 2. 1 .3, the power corrections in g [equation (8.6)] are connected with the corrections of order 
to the perturbative coefficients E^j^ (see also appendix D.l). 

It is therefore interesting to numerically determine the corrections to the leading factorial growth of the perturbative 
coefficients. We start with the available results for the first 300 perturbative coefficients and subtract the known 
leading asymptotics of the form Eq ^ —3^ +^ K\/tt, as well as divide the results of the subtraction by the 
leading asymptotics. An approximation for the coefficient of the correction term of relative order 1/K can be 
found via multiplication of the results of the previous operation by K. This approximation can be improved by 
elimination of power corrections of higher order in 1/K via "extrapolation to K ^ oo." For this latter step, 
various algorithms may be used. One possibility consists in the Neville algorithm [14], which is a variant of 
Richardson extrapolation [15]. The result is again a numerical estimate for the coefficient of the 1/K-Xenn. The 
exact, rational form can then be found easily by routines built into modern-day computer algebra systems, for 
example the Rationalize function of [16]. These operations may be repeated, and conjectures may be found for the 
coefficients multiplying the corrections of relative order 1 /K to the factorial growth of the perturbative coefficients. 
We find 



lK+1 



Kl 



1 



53 1 
18 K 



1277 1 
648 K^ 



405 395 1 
34 992 K^ 



218 793 923 1 
2 519424 K^ 



35 929 260 709 J_ 
45 349 632 



O 



1 
1^ 



(8.7) 



If we rewrite the corrections to the leading factorial growth of the perturbative coefficients in terms of a factorial 
series. 



77.0 



3^+1 K\ 

TT 

3^+1 K\ 

IT 



1 - F 7 



1-^ 



02 



K K{K-1) 



as 



(8.8) 



K{K-l){K-2) K{K- 1) {K-2){K-2,) " 

then it is easy to identify the coefficients entering into (8.6) with the Uj coefficients. Based on the discussion in 
chapter 2.1.3, it is easy to show that 

3^' (-a,) = eo,2y . (8.9) 
The results listed in equation (8.7) therefore lead to the following conjectures for the two-instanton coefficients: 



60,211 — — 



60,214 = — - 



53 



6 ' 
141 158 555 

31104 



60,212 — — - 



1277 



72 



60,213 — ~- 



336437 



1296 



60,215 = — - 



17 542 610 737 
186 624 ' 



(8.10) 
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The numbers 53, 1277 and 336 437 are prime. 



8.3 Perturbative Expansion and Instanton Theory 

We will now use a completely different and more direct route to the calculation of instanton coefficients: It is based 
on the explicit solution of the quantization condition (2.21) by an ansatz as given by the resurgent expansion (2.13). 
In order to carry out a calculation, we first have to determine the coefficients in the expansions of B^^{E, g) and 
^dw {E, g) defined in (2.23a) and (2.23b). An accurate calculation requires more terms than those given in (2.24a) 
and (2.24b). 

The calculation of the perturbative function Bi^{E, g) can be carried out easily based on equation (3.31). Indeed, 
recursive algorithms are known [8]. The generalization of (2.24a) to higher orders is as follows: 



Sdw(B,s) = E + g 



V 2 4 32 J 

. / 45045 5 45045 3 31185 
+g -E H E-^ H 

V 4 8 64 

5 / 969969 B 2807805 4 1924923 2 159159 \ 
+9 + E* H E^ H + 

V 4 16 64 256 ) 

6 / 22309287 7 88267179 5 100553453 3 25746721 \ 
+9 H E^ H E^ H E | 

" V 4 16 64 256 y 

7 / 2151252675 g 2788660875 „ 9526065549 4 
+0^ H E'' H E* 

V 16 16 128 

2506538463 2 692049787\ 

H H 

256 4096 / 



8 / 214886239425 „ 353522522925 7 1691601686775 5 

+g B H E^ H E^ 

V 64 64 512 

7583829646255 , 663834081625 



1024 16384 



E] +0(/). (8.11) 



The calculation of the "instanton function" A^-^ (i?, (/) is a httle more difficult (see the appendices F.6 and F.7): The 
calculation is based on the evaluation of successively higher orders of the contour integral of the WKB expansion 
as given in (3.57). The generalization of (2.24b) to higher orders is given by: 



3g \ 2 

3 187 ■ 
227 E"^ + — E 

^ ,'47431 4 34121 ^ 
+9 —TTT- E + E 



12 24 576 

, ,'317629 „5 264725 ^3 842909 \ 

+g E^ + E^ + E + 

'4 6 192 J 

26145967 « 16601579 4 63996919 2 6167719 \ 

E'' H E* H E^ H 

15 12 240 960 J 

e / 812725953 7 3490889111 5 
+g E^ + E^ 

V 20 80 

4398906487 3 1280980929 

H E'^ H E 

320 1280 

7 / 443323117271 » 265222473925 « 4948336000477 4 
E" H E" H E 

V 448 192 7680 

10166658134543 , 3228992367577 \ 

H E^ H + 

107520 172032 / 

o / 22315986340103 9 4909541135621 7 29042282605297 5 

+g + E^ H E^ 

" V 890 112 1024 

15025619362858 3 29626122220400 \ g, 

H E H E + O(g^) . (8.12) 

21504 688128 y t va / 
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8.4 Instanton Coefficients 



All coefficients, up to eight-instanton order, up to seventh order in g, have been calculated (see Table F.l). This 
requires the WKB expansion up to 5*8 (5 , -B , g) . Due to the pref actor 1 /g, we obtain relevant terms with positive 
powers of E from 6*10 of order . This is irrelevant if we want to calculate all terms up to the order of g^, and we 

may therefore neglect 6*10 (see also appendix F.6). 

The perturbation theory function Bdv,{E,g), together with the instanton function Advf{E,g), entirely determine 
the perturbative expansion about the n-instanton effect (n arbitrary) to order g^, i.e. all coefficients eo,nki with 
I < 8 in the resurgent expansion (2.13), for arbitrary n and k. In order to illustrate the effect of the instanton 
contributions on the energy levels, we give here formulas for the complete instanton effects through the order of g^ 
rather than fists of coefficients. This is inspired by the notation (2.15). For the states with parity e = ± and N = 0, 
the one-instanton effect reads [see also (2.27) and (2.28)]: 



= (l-T^.9-=.9^ 



71 6299 o 2691107 



12" 288 " 10368 
2125346615 4 509978166739 846134074443319 



9 9 A^r^r^c.^r.,^ 9 



497664 " 5971968 ^ 429981696 

_ 262352915454007369 7 _ 717976540715437267525 g 

5159780352 ^ 495338913792 ^ 

+0(/)) . (8.13a) 

Note that according to [8], many more terms in the perturbative expansion about one instanton may be calculated. 
There are even recursive formulae available. The first coefficients in the perturbative expansion about one instanton 
therefore read as follows: 

71 6299 
eo.ioo = 1, eo,ioi = -Y^, eo,io2 = -^^, 

2691107 2125346615 
eo,io3 - -^0368"' " 497664 " ^^'^^^^ 

The coefficients eo,ioi {1 = 0,..., 300) are available in electronic form [7]. 

Now we consider the two-instanton energy shift. Equation (8.5) clarifies that the two-instanton shift is given by 
the sum of two infinite series in g, which differ by the presence (or absence) of the logarithmic prefactor xig)- The 
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two-instanton effect is parity independent, and the explicit formula for A?^ = reads 



eig) 



, . 53 1277 2 

X{9) ( 1 - y - 9 - 



■9 



141158555 

31104 
14922996684685 g 

6718464 ^ 
3279840218627988925 
1934917632 
23 53 

y ~ "6 

45941 336437 
■7 



72 

17542610737 



336437 
1296 ■ 



186624 " 

2359159111315567 7 

40310784 ^ 



9' + 0{9' 



+ 7+ -^-ir7 



■13 
+ '12 



1277 



7 9 



72 

20772221 



144 1296 
12783531515 17542610737 



+ - 



+ - 



+ 



62208 
2110670726275 

373248 
2249987155449805 

13436928 
429051019455941467 



2592 
5 



141158555 
31104 



7 9 



7 9 



186624 

14922996684685 



7 9 



6 



6718464 
2359159111315567 



7 9 



80621568 



40310784 

3279840218627988925 
1934917632 



7 5 + 



(8.14a) 



The logarithmic coefficients 60,21; il = 0, . . . , 8), determined here by direct analytic calculation, verify the con- 
jectures {I = 0, . . . , 5) which were previously derived based on the corrections to the large-order growth of the 
perturbative coefficients (8.10) and the dispersion relations discussed in chapter 2.1.3. 

The three-instanton correction to energy eigenvalues, for states with N = 0, involves three powers of the nonper- 
turbative factor ^(g). The analytic expressions become very involved, as there are three infinite perturbative series 
in g, multiplying the logarithmic terms x^(5)' xCff)' ^^id a series multiplying the terms which lack the logarithm. 
The coefficients of the instanton expansion can still be calculated in closed analytic form, by direct reference to 
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the quantization condition (2.22). Up to terms of order g^, the three-instanton shift is given by 

£;(3)/ ^ ^ _ c3f ) [ 2(- ) / i _ Mi _ 1?!^ ^ _ 89635 3 _ 27325159 4 _ 2362765483 5 
€,0(9) — « 6 (.Sj (.5) I 2 8 ^ 64 ® 256 ^ 4096 ^ 16384 ^ 

1378627835215 « 148426925165305 7 139916898211047653 » 9,1 

s g g + 0(g ) } 

393216 1572864 50331648 J 

, , f / 63 141 \ /825 489 \ 2 / 24483 89635 \ 3 

+x(9)|37+(---— 7J -— 7J S +(-—-— 7) / 



/ 3369081 27325159 \ 4 , f 1534609037 2362765483 5 



256 2048 



, 179089589633 1378627835215 . , 
' 16384 196608 ' ^ 



\ 4 / 1534609037 2362765483 \ 
'^J ^ \ 4096 8192 V 



44105311958423 148426925165305 

131072 786432 / 

17311037863755515 139916898211047653 \ g 

7 9 +0(g )\ 

1572864 25165824 / ^ 'J 

f / 3 2 1 . A / 17 63 141 2 47 , ,^ 



2 

9 



/ 297 825 489 2 163 , , 
+ 1 7 7 C 2 

V 2 8 64 64 ^ 

/ 21199 24483 89635 2 89635 , , 
+ 7 7 C(2) 

V 192 64 256 768 

/ 271647 3369081 27325159 2 27325159 \ 4 

128 256 4096 12288 ^ ') ^ 

( 9180880297 1534609037 2362765483 2 2362765483 5 
+ 7 7 C(2) 

V 61440 4096 16384 49152 ' ) " 

( 572797894871 179089589633 1378627835215 n 1378627835215 , \ e 
+ 7 7 C 2 g^ 

V 92160 16384 393216 1179648 ''^ ' J " 

9833533558034329 44105311958423 

41287680 131072 

148426925165305 2 148426925165305 , A ^ 

7 C(2) g 

1572864 4718592 ' j " 

/ 150070370854927511 17311037863755515 

V 16515072 1572864 

139916898211047653 2 139916898211047653 , ,\ « g^ll 

7 C(2) S +0(s )H ■ (8-15) 

50331648 150994944 •'^ ' J " ^ ' }\ 



8.5 Reference Values 

All accurate verifications of the theory of instantons, at small coupUng g, necessitate a very precise determination 
of energy eigenvalues in the small-g region. A numerically accurate determination of eigenvalues can be a rather 
hard problem, even in a one-dimensional case. However, calculations are simplified when using multiprecision 
libraries [17-19]. In the current chapter, we intend to give reference values for specific small-;; example cases. 
The results for -£±.0 in all cases are rather close together, and the tiny energy differences then point toward the 
instanton-mediated energy shifts. In [9], we indicated two such values, valid up to 180 decimals, for the case 
g = 0.001. Here, we restrict the discussion to slightly fewer decimals while stressing that essentially arbitrary 
accuracy is available when using appropriate numerical algorithms [9]. For g = 0.002, we obtain for the ground 
state to 100 decimals 

i;+,o (0.002) = 0. 49798 16336 05614 52785 33444 97756 93929 30135 47830 
41905 18141 65406 30388 85981 28620 52208 12662 00662 , 



(8.16) 



whereas the first excited state has the energy 



£;_,o(0.002) = 0. 49798 16336 05614 52785 33444 97756 93930 90653 58949 
74478 49607 06416 37435 43472 00173 52993 86232 45517 . 



(8.17) 



The two energies differ on the level of the one-instanton contribution which is given by equation (8.13) and has 
opposite sign for states with opposite parity. Indeed, the energy difference i?_,o(0.002) — iJ+,o(0.002) can be cal- 
culated to high accuracy by simply evaluating the partial sums of the perturbative expansion about one instanton, in 
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analogy to the case g = 0.001 considered in [9]. As mentioned previously, the coefficients eo,ioi (/ = 0, . . . , 300) 
are available electronically [7] to the interested reader. For g = 0.005, we obtain as reference values, to 45 deci- 
mals, 

i;+,o(0.005) = 0.49488 15073 206 47 62721 54849 56033 00402 87323 00532 , (8.18a) 

and 

i;_,o(0.005) = 0.49488 15073 206 99 29084 80981 89004 43542 33907 22906 . (8.18b) 
The mean energy is 

j;+,o(0.005) + j;-,o(0.005) ^ p 49488 15073 20573 45903 17915 725I8 . (8.18c) 

Again, the real part of the Borel sum of the perturbation series (2.50) is slightly different, due to the two-instanton 
correction. 

Re B |£;^°Vo.005'l | = 0.49488 15073 20673 45903 17915 68107 . (8.18d) 
The energy difference is 

i;_,o(0.005) - £;+,o(0.005) = 5.16636 32613 23297 14313 9465 8 42237 x IQ-^'' . (S.lSe) 

Due to the three-instanton shift, the quantity i?_ o(0.005) — £^+^0(0. 005) is not quite equal to (twice) the real part 
of the Borel sum of the perturbative series (2.51) about one-instanton, 

2B^E^^\om5)^ = 5.16636 32613 23297 14313 9465 6 65773 x 10"". (8.18f) 

The case g = 0.007 yields the following energies 

£;+,o (0.007) = 0.49276 251 38 34552 88807 85274 35274 12515 77733 55043 , (8.19a) 

and 

i;_,o(0.007) = 0.49276 251 44 24291 38099 89370 58280 41759 92080 41287 . (8.19b) 
The mean energy is 

£;+,o(0.007) + g-,o(0.007) ^ 0.49276 2514129422 13 453 87322 . (8.19c) 

The real part of the Borel sum of the perturbation series (2.50) is sUghtly different, due to the two-instanton 
correction. 

Re ^1^;^°^ (0.007) I = 0.49276 25141 29422 13 399 23438 . (8.19d) 

The energy difference is 

£;_,o(0.007) - i;+,o(0.007) = 5.89738 49292 04096 2 3006 29244 14346 x 10"^° . (8.19e) 

Again, due to the three-instanton shift, the quantity E^^ (0.007) — i?+,o (0.007) is slightly different from the real 
part of twice the Borel sum of the perturbative series (2.51) about one-instanton, 

2 (0.007) I = 5.89738 49292 04096 2 0714 40266 74322 x 10"^°. (8.19f) 

Finally, for the slightly less problematic case of g = 0.01, we only list the energies 

i;+,o(0.01) = 0.489497520 976030, i;_,i(0.01) =0.489498132 721197. (8.20) 
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Function A(g) for the Double-Well 



1.02 



0.0 



Numerical Data 

Asymptotic Expansion up to 0(g ) 




0.01 0.02 

g 



0.03 



Figure 8.1: Double- well potential: Comparison of numerical data obtained for the 
function A{g) with the sum of the terms of order (where j < 10) of its asymptotic 
expansion for g small, where we express both the numerator as well as the denomina- 
tor of (2.52) as a power series in g. [When the resulting fraction is further expanded 
in g, then the leading terms of order g and g^ are given in equation (2.53) and rep- 
resented graphically in figure 2.2]. The higher-order terms employed here are easily 
derived using equations (8.13) and (8.14). There is good agreement between numer- 
ically determined ("exact") values (data points) and the smooth curve given by the 
analytic asymptotics. 



8.6 The Function A{g) for States with N = 

We recall the definition of the function A{g), given in (2.52), 

{i {E+,, + E^,o)-B{E^"\g)]] 



A(.g) 



(£;+,o-S-,o) [ln(2.9-i)+7] 



(8.21) 



We provide here, based in part on the data given in chapter 8.5, some reference values of the for the function A{g), 
for a somewhat larger coupling g as compared to those listed in Table 2.1: 

A(O.Oll) = 1.01147, A(0.012) = 1.01211, A(0.014) = 1.01322, A(0.015) = 1.01368, 

A(0.016) = 1.01408, A(0.017) = 1.01439, A(0.018) = 1.01465, A(0.019) = 1.01482, 

A(0.020) = 1.01489, A(0.021) = 1.01486, A(0.022) = 1.01464, A(0.023) = 1.01444, 

A(0.024) = 1.01392, A(0.025) = 1.01339, A(0.026) = 1.01233, A(0.027) = 1.01144, 

A(0.028) = 1.00991, A(0.029) = 1.00816, A(0.030) = 1.00653. (8.22) 

The available higher-order corrections (in g) to the two-instanton energy shift [see equation (8.14)] allow for 
an accurate comparison of the numerically determined energy eigenvalues with the asymptotic expansion of the 
function A(g) (see figure 8.1). We recall that the calculation of the instanton coefficients (8.14) relies on higher- 
order WKB expansions are described in appendices F.6 and F.7. 
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8.7 Leading Instanton Effects for States with N = 

Based on the simphfied quantization condition (4.60), derived using the path-integral formalism, it is relatively easy 
to determine the leading (in g) coefficients of the instanton expansion (2.13) to any order in the nonperturbative 
factor ^(9) = cxp( — 1/6(7) / V^fff^ zeroth order in g (this is equivalent to the coefficients (iN.nki with I = 0). In 
the current chapter, we intend to provide results for the states with = 0, up to eight-instanton order, in closed 
analytic form. The coefficients display an interesting analytic structure; a rapid growth of their absolute magnitude 
is observed in higher instanton orders for constant (zeroth) order in g (see also Table F. 1). We briefly recall here that 
the leading two- and three-instanton coefficients have been given in (2.40), (8.14) and (8.15). The four-instanton 
correction is discussed in appendix F.8. The four-instanton coefficients determine the polynomial P4 implicitly 
defined by (2.39): 

g 

eo,430 = g, 60,420 = 87, 60,410 = 87^ + 2^(2), 

eo,40o = ^7' + 27C(2) + ic(3). (8.23a) 

The five-instanton coefficients which enter into P5 , read: 

125 

60,540 = (8.24a) 
125 

60,530 = -^7, (8.24b) 
6 

125 ^ 25 , , 

60,520 = — 7' + yC(2), (8.24c) 

125 25 5 
60,510 = -^7' + y7C(2) + 3C(3), (8.24d) 

125 25 5 29 

60,500 = ^7' + y7'C(2)+3 7C(3) + y^C(4). (8.24e) 

The result for the six-instanton coefficients are a little more complex: 

54 

60,650 = -T- > (8.25a) 
5 

60,640 = 547, (8.25b) 

60,630 = 108 72 + 18C(2), (8.25c) 

60,620 = 108 73 + 547C(2) + 6C(3), (8.25d) 

60,610 = 5474-h5472C(2) + 12 7C(3) + ^C(4), (8.25e) 
54 

60,600 = y7' + 187'C(2) + 672C(3) 

+C(2)C(3) + §7C(4) + ^C(5). (8.25f) 
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The seven-instanton coefficients, to lowest order in g, read: 



720 

120 



(8.26a) 

7 , (8.26b) 



7° 2 2401 , , 

— 7 H C(2) , (8.26c) 

48 48 ^ 

7^ , 2401 , , 343 , , 

^7=^ + ^7a2) + ^f(3), (8.26d) 

7^ 4 2401 2 . 343 , , 1911 , , 

^7^ + — 7'C(2) + — 7C(3) + — C(4), (8.26e) 

^7^.^73C(2).515,= C(3).?C(2)C(3) 
1911 7 

+ -Tr7C(4) + -C(5), (8.26f) 
16 5 

6 , 2401 4 343 3 49 

— 7 +—7 C(2) + — 7 C(3) + -7C(2)C(3) 

7 1 Q 1 1 7 7SQ 

+ ^ C'(3) + 7' C(4) + - 7C(5) + — C(6) . (8.26g) 



The eight-instanton effect features the following coefficients: 



(8.27a) 

7 , (8.27b) 

2^^ 

7' + — C(2) , (8.27c) 



2" 

eo,870 = — 

2" 
45 
2" 

eo,850 = 7^ . ■ 15 

214 2" 

eo,840 = — 7^ + — 7C(2) + yC(3), (8.27d) 

704 

eo,830 = — 7'' + — 7''C(2) + — 7C(3) + — C(4), (8.27e) 

2" 5 2i= 3 2I0 2 2^ 

eo,820 = ^7' + — 7'C(2) + — 7 C(3) + yC(2)C(3) 

2= 

+704 7C(4) + — C(5), (8.27f) 

5 

eo,8io = ^7' + ^7*C(2) + ^7'C(3) + ^7a2)C(3) + ^C'(3) 

+704 7= C(4) + — 7C(5) + 62 C(6) , (8.27g) 

5 

eo.800 = Itt 7^ + ^ 7' C(2) + ^ 7' C(3) + ^ 7' C(2) C(3) 
315 15 9 3 

2^ 704 22 2^ 

+ V 7 C'(3) + — 7' C(4) + - C(3) C(4) + - 7' C(5) 
9 3 3 5 

+ 7C(2)C(5) + 62 7C(6) + ic(7). (8.27h) 
5 7 

It is intriguing to observe that in all individual terms contributing to a particular coefficient, the sum of the power 
of 7 and of all (integer) arguments of ( functions equals a constant. An analogous pattern may be observed in 
the context of specific sums that enter into the evaluation of higher-order corrections to the vacuum-polarization 
charge density around a nucleus, as used in the derivation of equations (66) and (67) of [20] (relevant formulas are 
given in appendix IV ibid.). 



8.8 The Function Ai{g) for States with iV = 1 

In this chapter, we will be concerned with the generahzation of the function A{g) [equation (2.52)] to states with 
= 1. The perturbation series for iV = 1, up to order g™, has been determined analytically based on the recursive 
techniques outlined in [8]. Coefficients are made available at [7]. The instanton coefficients for states with = 1 
are different from those hsted in equations (2.18), (2.19), (8.13), (8.13b) and (8.14), and in the Table F.l, where 
the states with AT = are considered. However, it is relatively easy to calculate the coefficients ei^nki with n <2 
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and / < 6 based on the quantization condition (2.22) and the explicit formulas for the function Bdw{E,g) and 
Adv/iE, g) given in (8.11) and (8.12). For the states with parity e and N = 1, the one-instanton effect is enhanced 
by a factor g~^ in comparison to the N — 0-states [see also (2.27) and (2.28)]: 



347 5317 o 15995159 ^ 28062012119 



12 288 10368 497664 " 
12918255230839 . 38332497252543415 . 7A 
^7196^ ^ 429981696 ^ + ""^^ ^ ) " ^'-''^ 

The coefficient ei_ioo is governed by the general formula (2.35). The two-instanton energy shift, for states with 
N = 1, reads 

(2) 4 2, , r , , A, 293 , 39823 2 2767481 3 902381531 4 

E, 1 = -rr £ (g) x(g) 1 9 H 9 9 9 

«,i g2 y g « 72 « 1296 " 31104 

233318679457 ^ 377478189958993 g , , 

9 9 + 0(g^ 

186624 6718464 

/ 61 293 \ / 21179 39823 \ 2 / 1592699 2767481 

^''-^^ + \T- — V '+\^ + ^V ' +1 — Sii li^^. 



548736791 902381531 \ i ( 33161353697 233318679457 

31104 31104 '< j 9 + y 93312 186624 

6148524786991 377478189958993 



9 



. (8.29) 



6718464 6718464 

Again, the coefficients 61^210 ™d 61^200 governed by the general formulas (2.37) and (2.38). 




Figure 8.2: Double-well potential: Comparison of numerical data obtained for the 
function Ai(g) defined in (8.35) with the sum of the terms up to the order of 
of its asymptotic expansion for g small, where we express both the numerator as well 
as the denominator of (8.35) as a power series in g. There is good agreement between 
numerically determined ("exact") values (data points) and the smooth curve given by 
the analytic asymptotics. See also Table 8.1. 



We consider the example g = 0.005, both in order to illustrate the calculation and in order to provide reference 
values for independent verification. The energy of the state (+, 1) as determined numerically reads: 

S+a (0.005) = 1.4632133515 77109 24646175934042775209 25202. (8.30) 
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The energy of the state (— , 1) is sUghtly higher: 

£;_,i(0.005) = 1.4632133515 9534145149 6975727023 6354918578. (8.31) 
The energy difference 

AE(0.005) = £;_,i(0.005) - E+,i(0.005) w 1.823 x 10"" . (8.32) 

The mean energy is 

(0.005) + (0.005) ^ i 45321 33515 86225 34897 93675 33725 . (8.33) 

The real part of the Borel sum of the perturbation series is 

B \^e[°\0.005)^ = 1.46321 33515 86225 34897 88950 50((1) . (8.34) 

Those decimal figures which agree are underhned. The difference is due to the two-instanton energy shift. The 
definition (2.52) of the function A which relates the one- and two-instanton effects to the Borel sum of the pertur- 
bation series, should be modified in order to accommodate for the changed structure of the leading two-instanton 
coefficients. We define the function Ai{g) as 

Ai((?) = 4-!^ 2 5^ ^. (8.35) 

(£;+,! -^_,i)'[ln(2ff-i)- 1+7] 

In order to determine the leading asymptotics of the function Ai{g), we have used the perturbative expansions 
about one- and two-instantons given in equations (8.28) and (8.29). If we additionally perform an expansion in 
inverse powers of the logarithm ln{g), we arrive at the result 

which generalizes the formula (2.53) to the case N = l.Of course, many more terms in the asymptotic expansion 
of Ai (g) can be determined based on equations (8.28) and (8.29). Sample values for the function Ai (g) are given 
in Table 8.1, and the function is plotted against its small-^ asymptotic expansion in figure 8.2. 



Table 8.1: The ratio Ai(g) as a function of g. The values obtained from the asymptotic ex- 
pansion for g small (up to the order g^", denoted "asymp.") are compared to numerically 
determined values (denoted as "num."). Useful numerical data for an independent verifica- 
tion is given in equations (8.18) and (8.19). 



coupling fi 


0.005 


O.OOG 


0.007 


0.008 


0.009 


0.010 


Ai{g) num. 
Ai{g) asymp. 


1.02097 
1.02098 


1.02434 
1.02434 


1.02748 
1.02748 


1.03039 
1.03039 


1.03304 
1.03305 


1.03540 
1.03542 



8.9 Three-Instanton Effects for AT = 

According to the equations (2.27), (2.39) and (2.40), the three-instanton shift of the states with = 0, at leading 
order in g, is given by 

Ei'\g) = -eeig) {P^M-V9)] + 0{gln\g))} , (8.37) 
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Figure 8.3: Double- well potential: Comparison of numerical data obtained for the 
function D{g) defined in (8.42) with the sum of the terms up to the order of 
of its asymptotic expansion for g small. Again, we express both the numerator as well 
as the denominator of the expression defining D{g) [see equation (8.42)] as a power 
series in g. 



where the expression [ln(— 2/5)] reads 



In - 



1 2 



+ 7 



TT 

12 



(8.38) 



In order to evaluate the physically relevant part of the three-instanton shift, we have to give an interpretation to the 
squared imaginary part generated by the logarithm, for positive g, under the analytic continuation of the logarithm 
to negative argument, 

In ^-^^ ^ In 0^ +i7r. (8.39) 

There is a potential problem. The explicit imaginary part of the three-instanton shift, generated by the analytic 
continuation of the logarithm, has to cancel the imaginary part of the (generalized) Borel sum of the one-instanton 
shift. In chapter 2.L3, this fact has been used in order to derive the leading large-order asymptotics of the pertur- 
bative expansion about one instanton. The imaginary part therefore cancels when the resurgent expansion (2. 13) is 
summed, and we need not consider it when evaluating real energy shifts. On squaring the right-hand side of (8.39), 
we isolate the physically relevant part of P^{g) which for g > reads: 



5>0 



^eP^{g) = - 









+7 







17 2 

TT 

12 



(8.40) 



According to the conjecture (2.13), the energy difference E_ fl{g) — E^ Q{g), for 5 > 0, is approximately equal 
to twice the value of the (positive) one-instanton energy shift of the state (—,0). The difference of the quantity 
E-fi{g) — E+fi{g) and the real part of the Borel sum of the perturbative expansion about one instanton is approx- 
imately given by the three-instanton shift. 



E_^,{9)-E+^o{9)-2B[E[°\g)} 



2E%{g) 



(8.41) 
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8.9 Three-Instanton Effects for AT = 



Table 8.2: The ratio D{g) as a function of g. The values obtained from the asymptotic ex- 
pansion for g small (up to the order g^, denoted "asymp.") are compared to numerically 
determined values (see chapter 8.5, values are denoted as "num."). The interested reader may 
verify the numerical values based on the data given in equations (8.18) and (8.19). 



coupling g 


0.005 


0.006 


0.007 


0.008 


0.009 


0.010 


D{g) num. 
D{g) asymp. 


1.00883 
1.00877 


1.00931 
1.00920 


1.00944 
1.00925 


1.00919 
1.00891 


1.00856 
1.00816 


1.00754 
1.00698 



where we neglect the five-instanton contribution and recall that E_ q (5) > for (/ > 0. In order to define a function 
D{g) in distant analogy to (2.52) and (8.35), we should normalize to unity in the limit 5 — > 0. This can be achieved 
as follows: According to (2.35), the quantity i?_ o(.9) — E+,Q{g) is approximately equal to 2^(g) in the small-g 
limit, where ^((?) is defined in (2.27). The third power of E-fl{g) — E+fl{g) then compensates the three powers 
of ^{g) comprised by the three-instanton effect (8.15). In order to normalize to unity, the only remaining missing 
terms are a prefactor 4 and an inverse normalization factor ReP3 (g') which is given in (8.40). We then define the 
function D{g) as 

4 E_,o-E+,a-2B[E[°\g)} 



D{9) 



In ( - ) + 7 



17 o iE.,o-E+,oy 
12 



(8.42) 



.5. 

The leading asymptotics for small g read [we also expand in inverse powers of the logarithm, as in (2.53) 
and (8.36)]: 

Of course, more terms in the asymptotic expansion of D{g) for small g can be calculated in a straightforward 
manner on the basis of equations (8.13) and (8.15). A table of numerically calculated values of D{g), and of the 
values suggested by the expansion up to the order g^, is given in table 8.2, and a graphical representation is obtained 
in figure 8.3. In the obtaining the table 8.2 and figure 8.3, the expansion of D{g) is carried out through the order 
of g^, and no additional expansion in inverse powers of the logarithm ln(2 g~^) is performed. 



36 



Chapter 9 

Conclusions 



The central theme of the current article is the discussion of (multi-)instanton effects, which manifest themselves in 
nonperturbative, nonanalytic (in the coupling) contributions to the path integral which in turn defines the partition 
function (see chapters 1, 2 and 3). The expansion about nontrivial saddle points of the Euclidean action leads 
naturally to energy shifts which involve nonanalytic factors of the form exp(— a/g). The instanton interaction and 
multi-instanton effects find a natural interpretation in terms of classical trajectories along which the particle may 
oscillate between degenerate minima. A particular detailed discussion of instantons in the (symmetric) double-well 
problem is provided in chapter 4, and the considerations are generahzed to asymmetric wells in chapter 5, to the 
periodic cosine potential in chapter 6, to radially symmetric oscillators (self-adjoint extension and resonances), 
and potentials of the Fokker-Planck type (chapter 7). In chapter 8, we complement the preceding chapters by 
higher-order analytic as well as numerical calculations related to the double-well problem, up to the order of eight 
instantons (see also appendix F). 

The conjectured quantization conditions are presented in chapter 2.3. For completeness, we give a Ust of the classes 
of potentials for which conjectures on quantization conditions are discussed in the current article, 

• the double- well potential (section 2.1 and chapter 4), 

• more general synnmetric potentials with degenerate minima (section 2.2.1 and appendix E), 

• a potential with two equal minima but asynnmetric wells (section 2.2.2 and chapter 5), 

• a periodic-cosine potential (section 2.2.3 and chapter 6), 

• resonances of the C'(i^)-symmetric anharmonic oscillator, for negative couphng (sections 2.2.4 and 7.2), 

• a special potential which has the property that the perturbative expansion of the ground-state energy vanishes 
to all orders of the couphng constant (see chapters 2.2.4 and 7.2.3), 

• and eigenvalues of the 0(!/)-synunetric anharmonic oscillator, for negative couphng but with the Hamilto- 
nian endowed with nonstandard boundary conditions (section 2.2.5 and 7.3). 

Indeed, exact results may be derived for large classes of analytic potentials, as discussed in chapter 2. The consider- 
ations presented in chapters 3 .2 and 3 .3 are important for more general investigations in later chapters of the current 
work, and indeed might be useful for any further conceivable generahzations. Specifically, the rather general quan- 
tization condition summarized in equations (3.55) and (3.57) may be adapted to large classes of analytic potentials. 
The perturbative B-function and the instanton j4-function are known to higher orders in the couphng g, for a many 
of the potentials discussed here. For the periodic cosine potential [Eqs. (6.31) and (6.32)] and the Fokker-Planck 
potential [Eqs. (7.32) and (7.35)], as well as the 0{v) anharmonic oscillator [Eqs. (7.48) and (7.49)], we perform 
higher-order calculations of the B- and A-functions up to the order g'^. For the double-weU potential, we perform 
calculations up to 0{g^) [see Eqs. (8.11) and (8.12)], i.e. up to eight-instanton order. 

The investigation of nontrivial saddle points of the Euclidean action allows for the identification of finite-action 
contributions to the path integral, the instantons, which may be calculated exactly in the hmit of a large separation 
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of the instanton tunneling configurations (see the discussion in chapter 4.4.2). Approximate quantization conditions 
may be derived on the basis of related considerations that evaluate the sum of the leading instanton effects, with an 
arbitrary number of oscillations on the instanton trajectory but considering only the leading term in the instanton 
interaction [see equations (4.60), (5.16), (6.29), (7.22), (7.29), (7.66), and (7.68)]. 

We now briefly recaU general consequences of the resurgent expansion for the energy levels of the double-weU 
potential (2.13), which reads 



It is thus clear that ordinary perturbation theory, which reUes on an expansion in powers of the coupUng constant, 
is not even qualitatively sufficient for a description of the energy levels as a function of the coupling constant g. 
In other words, whenever there are nontrivial saddle points of the classical action, it is not sufficient to consider 
perturbation theory about the trivial saddle point of the action. Generalized expansions are required which typically 
include nonanalytic factors of the form exp(— a/gf) where [see equation (2.60)] 

/■90 

= 2 / dg V2^(5) , (9.2) 

JO 

q and go are the positions of the degenerate minima, and V is the potential. 

Rather important differences prevail among cases where there is parity symmetry, and asymmetric cases. Indeed, 
for the (symmetric) double-well potential, the splitting of opposite-parity states is given by a one-instanton effect 
of the order of 2 exp[— a/(2g()]. By contrast, there is no such degeneracy of states in the case of asymmetric 
wells, and for the ground-state energy (7.31) of the Fokker-Planck Hamiltonian, we have an energy of the order 
of cxp(— a/g). For the asymmetric case, the classical trajectories have to return to the original minimum (this 
corresponds to an even number of "tunnelings" between the degenerate minima), or else the trajectories do not 
contribute to the path integral. This is discussed in chapters 4.2.2 and 5.1. 

Let us outline a few unsolved questions related to the discussed problems, which may be investigated in the future. 

• First and foremost, the exact WKB- and instanton-inspired methods as well as the numerical investigations 
discussed here may be extended to more general potentials [see the general equation (3.57)]. One example 
is provided by potentials with three or more degenerate minima which interpolate between the cases of two 
degenerate minima on the one hand and the periodic cosine potential on the other hand. A further, particularly 
intriguing case is the Fokker-Planck potential (see chapter 7.2.3) for which the perturbation series vanishes 
to all orders (the manifestly positive ground-state energy is determined in this case by instanton effects). 



Second, there is a certain unsolved issue related to the complete (re-)summation of the instanton expansion. 
Indeed, the nonperturbative factor exp(— a/5) may assume rather large values for moderate values of g. In 
the specific case (2.13), while the power series 

00 

^eN,nkid (9.3) 

(=0 

find a natural summation procedure in terms of the Borel method in complex directions, the large-order 
properties of the expansion in powers of 

(9.4) 



■ng 

are largely unknown. The explicit calculations in chapter 8 might be helpful for a verification of large-order 
estimates [regarding the expansion in exp(— l/6£f)]. The same applies to the entries of Table F.l. 

Third, as outlined in chapter 3.2.5, there is an interesting connection to fundamental properties of Borel 
summability. Note that even in situations where the perturbative expansion of eigenvalues is Borel summable, 
it is not clear whether the functions Bi(E, g) and A{E, g) [see equation (3.57)] are Borel summable in 51 at S 
fixed. Indeed, the wave function ^(g) is unambiguously defined only when the quantization condition (3.55) 
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is satisfied with N a non-negative integer. When E is not an eigenvalue, the solution of the Schrodinger 
equation is an undefined linear combination of two particular solutions. Thus, it may be interesting to in- 
vestigate the behaviour of the functions Bi{E, g) and A{E, g) when E approaches an eigenvalue, in closer 
detail. It is conceivable that this approach might reveal connections between the A- and the S-functions and 
eventually lead to a unified understanding of the structure of these functions. Again, the explicit calculations 
in chapter 8 might be useful in the context of a verification of such relations, if they indeed exist. 

• Fourth, there might be some possibilities for an extension of the exact methods discussed in the current 
article to higher-dimensional scenarios where our current understanding is more Umited. There is a well- 
known analogy between a one-dimensional field theory and one-dimensional quantum mechanics, the one- 
dimensional field configurations being associated with the classical trajectory of the particle. Therefore, 
model problems derived from field-theoretic models via a reduction of the dimensionality may find a nat- 
ural and exact treatment via the semi-classical methods discussed in the current article. Indeed, the loop 
expansion in field theory corresponds to the semi-classical expansion [21, chapter 6]. 

We conclude with a few remarks on conceivable applications to field theories, and in particular non-Abelian gauge 
field theories. These questions follow naturally after the famous examples found by Lipatov [22-24] for scalar 
theories. For a study of the non-Abelian case, one could be tempted to assume that the periodic potential discussed 
in chapter 6 might be a good starting point. However, one of the pre-eminent problems is related to the scale 
invariance of the classical equations of motion (classical chromodynamics) of the fields which are invariant under 
the transformation A^{x) — > A A^{\x). In a multi-dimensional field theory, instantons are relevant at all possible 
length scales (small and large separations), and it is not a priori clear how to perform the integration of the small- 
and large-scale instantons. An exphcit infrared cutoff destroys the gauge invariance of the theory. For a survey 
of field-theoretic aspects related to instantons, one may consult the following review articles [25-29] (the Ust is 
necessarily incomplete). The nonlinear a model (see e.g. [30, chapter 14]), or the (fA theory [30, chapter 34], or 
even manifestly nonrelativistic field theories might provide better defined scenarios for a systematic, desirably 
exact study of instanton effects. In general, a better and more systematic understanding of situations which have 
by now precluded a thorough analysis, might be gained by suitable generalizations of the methods discussed here. 
As demonstrated, these have been successfully applied to solvable model problems, in one-dimensional quantum 
mechanics, which is equivalent to a field theory (point-like in space, with one time dimension). 
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Appendix F 



WKB Calculation and Instanton 
Contributions 



F.l An Important Asymptotic Expansion 

To obtain the perturbative expansion of the function A{E,g) (equation (3.57)) we need the WKB expansion 
(B.13) of S{q). To extract A from equation (3.57) we have first to expand the term — lnr[i — B{E,g)] + 
B{E,g) \n{—g/2) for B large. This necessity arises because B{E,g) •-^ E = 0{g^^) in the context of the 
WKB expansion. 

For the derivative V'' of the V' function [^(2:) = d/dz In r(2;)], the following integral representation is known: 

^'(\ + ^= \ r \ e-*Mi = |1 r f i + zV (F.l) 

V2 ) 2 io sinh(i/2) dz^ \2 / 

For large z, the integral is dominated by the region of small t, and one may therefore derive an asymptotic expan- 
sion of + z) by expanding the integrand [excluding the exponential factor exp(— t z)] in t, and subsequent 
integration of the terms resulting from this expansion. Since {z) = d'^/dz'^ Inr(z), the integration constants 
have to be adjusted properly. One obtains the following asymptotic expansion of the F-function, 

/I \ 1 17 

Inr - + z = zlnz-z+-ln{2Tr)- — + 



2 J 2 ^ ' 2^z 2880^3 



40320^5 215040 608256^7 

F.2 Properties of the Mellin Transform 

The MeUin transform of a function F is commonly defined as 



M(s) 



I AEF{E)E-^-^ . (F.3) 
Jo 



A rather useful property of the MeUin transform is this: If the function F can be expanded into a power series in 
its argument, then the transform M develops poles at positive integer argument, and the corresponding residues 
yield the coefficients of the power series for F. Moreover, for functions F{E) which contain logarithmic terms 
of the form E^ In E, the MelUn transform has double poles whose residues, again, give the coefficients of the 
logarithmic terms. 

Let us consider a slight generalization of the Mellin transform with an arbitrary upper linnit of integration L. Then, 



Jo 



^ jn—s 1 

d£;£;"£;-^-i = = +lnL + 0(n-s). (F.4) 
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F.3 Mellin Transform of the WKB Expansion 



The residue at s = n is — a. For a term of the form 



L T ii-s T n-s 1 r 

dEE-lnEE-^-' = ^ L InL 







(n — s)2 n — s 
1 In^L 



The coefficient of the leading term in the Laurent series about n = s is independent of L. This is important for the 
considerations presented in the following chapters. 

F.3 Mellin Transform of the WKB Expansion 

To calculate the successive terms of A{E, g) in an expansion in powers of g we have then to replace in equation 
(3.57) S+ by its WKB expansion, and expand each term for Eg, that is E small. A standard method to obtain 
this expansion which contains only powers of E and powers of E multiplied by In E is to calculate the MelUn 
transform of integral (3.57). We thus consider the function 

M{s) = ^ jdEdqS+{q,E,g)E-'-'^ . (F.6) 

One verifies, replacing 5*+ by its WKB expansion, that the function M{s) has double poles at integer values of s. 
The residue of the double pole at s = n yields the coefficient of —E^ \nE in the expansion of the integral (3.57) 
for g small and the residue of the simple pole the coefficient of — 

^ E", (F.7a) 

n — s 

2 

ln£;. (F7b) 



1 



The formulas lead to a proper identification of the terms in the Mellin transform of the contour integrals of succes- 
sive orders of the WKB expansion, with the expansion in g of the contour integrals themselves of the successive 
orders of the WKB expansion. Therefore, the equation (F.7), for each order in the expansion in g, is in fact a Mellin 
backtransformation. 

We now investigate the contour integrals of successive orders of the perturbative expansion of the WKB expansion 
(sic!). We recall that the WKB expansion is an expansion in powers of g at gE fixed [equation (B.13)], and that a 
general result for So has been given in (3.40). [A general second-order result for 52 has been given in (B.16).] 
Expansions in g of 5*0 and 52 are given in (B.18) and (B.20), respectively. Formulas relevant for the contour 
integrals are given in (B.19) and (B.22). 

We call Io{s) the leading contribution to M (,s) coming from Sq, and /2(s) the contribution coming from g^ S2- 

M(s)=/oGs)+/2(s)+0(.g-'). (F.8) 
The function /o(s) is then given by [see equations (B.18) and (B.19)]: 



/o(s) = -^^j>dEdq[2gE-U^{q)f'^E-''-^ 

= _^ j^Edq [2gE-U\q)f^ (gE)-'-' 

where q = and q = go are the two minima of the potential. The function /o(s) has a simple pole for s = and 
double poles for s positive integer, indicating the presence of logarithms in higher powers of g. 
The contribution from g^ S2 is [see equation (B.21)]: 

= i(2»)'«r(-„pIfcV2_ jy,u-'-'-u-. (F.io, 

The function Ii{s) has a simple pole for s = — 1 and double poles for s integer, s > 0. 
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F.4 General Potentials at Leading WKB Order 



F.4 General Potentials at Leading WKB Order 



The residue at s = of /o(s) as defined in (F.9) is proportional to the instanton action (5.6): 



rQo 
Jo ' 



dqUiq) = 



(FIX) 



The residue and the double pole at s = 1 are also of crucial importance for the quantization condition (2.64). 
Indeed, by calculating the residues of the double and simple poles of Io{s) at s = 1, we obtain the coefficients 
of ElnE and E and thus terms which are generated only by the expansion of the F-functions and and the factor 
{-2C^/g)^ of equation (3.57). We now combine (B.18), (B.19), and (F.9) in order to evaluate /o(s). If we denote 
by 1 and the values of the second derivatives of the potential at the minima q = and q = qo, respectively [see 
(2.70)], then 



'/ ' 

Jo 



dqU'-^^{q) 



= 2 dq 
Jo 



U'-'^iq) 



^l-2s 



(w(go - q)y 



1 4- / il-2« 



(l-s) 



2 / dq 

Jo 

1 



1 



U-\q)-- 



1 



q u{qo - q) 



, ,1 + - +2 1 + - lngo + 2 . 

(1 — S) V UJ ) \ UU ) UJ 



(F.12) 



It follows, according to the correspondence (F.7), that the contributions to the integral on the left-hand side of 
(3.57) are 



E 



{\n{-E/w) -l)+E {\n{-E) - l) 



1 + - ) ln(-c?/2)-2— -21nC-2 (1 + - lingo 

iO I LO \ LO ' 



with 



90 

lnC= / dq 

10 



1 



1 



U{q) q w{qo-q) 



This result is fully consistent with equations (2.73) and (5.3) as well as (C.IO) below, 

B,{E,g)=E + 0{g), B^i^E^g) = E/uo + 0{g), ^ ql {u;^ C^^) 



i/(i+w) 



(F.13) 



(F14) 



(F15) 



We have used here the relation B{E,g) = E [1 + 1/lo) which is vahd for asymmetric wells in the convention 
(2.70). The result implied by this calculation, if combined with equations (3.55) and (3.57), is fully consistent with 
the quantization condition (5.16), which in turn is the expansion of the quantization condition(2.64) at leading 
order in g. 



F.5 Symmetric Potentials: Next-to-Leading WKB Order 



We now specialize the treatment of the previous chapter to synmietric potentials. In the case of symmetric potentials 
with degenerate minima, which have loi = LO2 = 1 in the sense of (2.63) and (2.70), one can introduce the 
parameterization (B.6). Then, 



dqU^-^\q) = 21 du 

Jo \J pyu) 

qo pu-o 

3-2sr„\Trr2r„\ _ o / J„,„,-3-2s 



Uo 



u 



l-2s 



dqU-^-'\q) U'\q) 



= 2 duu~ 
Jo 



(F16a) 
(F16b) 
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F.5 Symmetric Potentials: Next-to-Leading WKB Order 



where uq is the zero of the function p: 

p{uo) = . 

We now consider the residues of /q (s)- 

• The residue at s = of /q (s) yields the leading contribution to A{E, g) [equation (F.9)]: 



(F.17) 



a _ 2 
9 9 Jo 



dqU{q) = - 

9 Jo 



du 



(F.18) 



The residues of the poles at s = 1 of /o (s) require slightly more work: 



Jo 



du ■ 



1 



1 



VpR Jo u \^/p(uj 

We set, in accordance with (2.58), 

1 



1 +2lnuo + 0{s-l). 



lnC = 2 / — 



1 +21nuo. 



(F.19) 



(F20) 



lo u \y/p(u) 

Combining with the pole of the factor r(— s), one obtains the coefficients of E ln{E) and E, respectively, 

2ElnE-2E + 2Eln{g/2C). (F.21) 

We consider the residue at s = 2 of Jo(s). For the term of order g we need the small-u expansion (B.6) of 
the function p. Then, 



du ■ 



= 2 





+2 
1 



„l-2« _ ± Ql u2-23 

2 



n> du 

Jo 



1 



- 1 + ^a\u - I -ira2 + ^af 



with 



a2,2 



zOi2 + qCCi 



du 



1 ai 



1 

2-s 



a2,2 + 0{s - 2) 



(F.22) 



1 , 1 
—==-l + -aiu 

\'p{u) ^ 



1 



1 ^ 2 \ 2 

■-^a2 + gOl ) M 



H 1-2 — - ^2 + 77 a, Inuo . 

ug uo ' ° ' 



The factor then yields 



2«+i s-i r(-s)r(3/2) 



r(3/2-s) 2-s ' 2 

Combining factors, one obtains a contribution to gE^ In E: 



^ +l+g\n{g/2) + 0{s-2). 



1 3 
--a2 + -al]gE^\nE 



and to gE^: 



(F23) 



(F.24) 



(F.25) 



(F.26) 
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We now consider the residues of his). We first rewrite the expression (F.IO) in terms of the function (B.6) as 

/ duu 
Jo 



-3-2s 



(F27) 



The first residues of /2 (s) may now be evaluated as follows: 

• From s —1, one obtains the contribution — l/12i!^. 

• For s ^ 0, the factor in front of the integral in (F. 10) has the expansion 

_A_^ln(,/2)-l^. 

2.S 2 ' 3 

In combining (F.16b) with the parameterization (B.6), we conclude that the integral itself yields 



(F.28) 



/ 

Jo 



-3-2s 



du 'i 



-3-2s 



1 + - aiu - 



du 



= l--+lnno 



2«2-gai)- 



2ul 



r° du r /—- ^ 1 

Jo u I 2 



ai_ 
2 Wo 
1 1 



-ai u 



Using (F.28) and (F.7b), we infer that the contribution proportional to gin E thus is 

The coefficient of g (without logarithms) may be deduced from (F.28) and (F.7a), 

1 



(F29) 



(F.30) 



a2 — —Oil I u 



1 



ai 



4wo 4wo 



1 



1 



8 "^-32"^ 



ln(f)+^-21n.o 



(F.31) 



Meanwhile, we have evaluated certain contributions to the Melhn transform of the contour integrals of the leading 
orders of the WKB expansion. These Mellin transforms were denoted by Jq and I2, and after a MelUn back- 
transformation (F.7), these terms enter on the left-hand side of equation (3.57). We would now like to regroup 
these terms into a form which is amenable to the identification of the functions A{E,g) and B{E,g). These 
latter functions are found on the right-hand side of (3.57). Using the expansion of the F function and setting 
B{E, g) = E + gh2{E) + 0{g% one finds 



iln(27r)-lnrQ-i3(i?,5) 



B 



{E,g) In {-^\ 



2CJ 



(F32) 



Here, advantage has been taken of the fact that the WKB-expansion is an expansion in g at fixed Eg, wherefore 
E = 0{l/g). In order to obtain A{E, g), one has to subtract this contribution twice, because we have a symmetric 
potential, with two equal contributions from each of the two wells. One sees immediately that the terms of order 
5° cancel and as well as the term 2b2{E) \n{gE/2) [see equation (B.l 1)], as expected from general arguments. 
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We are now in the position to write down explicit expressions for the function A{E, g) at order g, for symmetric 
potentials which foUow the parameterization (B.6). The contribution at order g is then 



A{E, 9) = - +9 (^2,2 + a2,o) + 0{g^) , 



where 02,2 and 02,0 are given by 



^2,2 



0.2,0 



du 

2 Jo ^ 



--a2 + -a, 



02,2 , 



V p(u) - 1 - -aiu - 



1 1 2 \ 2 

- 0.2 - -«! ) M 



The expUcit result for 0,2,2 can be found in (F.23). 



(F33) 



(F34a) 



(F34b) 



F.6 Explicit evaluation for a Special Family of Potentials 



Let us remember at this stage that the basic paradigm for the evaluation of the "instanton function" A{E, g) that 
enters into (2.24b) is the following: start from the conjectured structure of the WKB expansion (3.57), 



1 



9 Jc 



AzS+{z) = A{E,g)+\n{2-K) 



Y,{\nT{\-Bi{E,g))+B,{E,g)\n{-g/2Ci)\ . 
i=l 



from which using (3.43) 



exp 



- I AzS+{z) 
9 Jc 



1 = 0, 



(F35) 



(F36) 



the quantization condition (2.64) foUows immediately. The perturbative expansions Bi{E,g) can in general be 
evaluated easily using the techniques described in chapter 3.2.5. The goal is the calculation of A{E, g). 

One may calculate successive orders in the WKB expansion using the algorithm (3.41), resulting in an approxima- 
tion for the left-hand side of (F.35). Because the Bi{E, g) on the right-hand side of (F.35) may easily be calculated 
and expanded according to (F.2), it is then possible to calculate A{E, g) by subtracting those terms that are gen- 
erated by the Bi{E, g) on the right-hand side of (F.35), from the result obtained for the WKB expansion on the 
left-hand side, and obtain a result for A{E, g). 

For classes of potentials for which the function p has a simple form, the expressions simplify. We thus consider 
now the class (B.23) of potentials which satisfy p{u) = 1 — 4u'". Then, the two integrals (F.16a) and (F.16b) can 
be calculated explicitly. One finds 



Jo 



^l~2s 

au = 

Vl - 4m" 



1 



,2(.-nM r(^(i-,s))r(i) 
r(i + ^(i-.))' 
%2(.+i)M r(-^(i + .))r(|) 
r(|-^(i + s)) 



The leading WKB order then yields 



1 



Us) = ±g^--L2'^i+i/m)+2-i/m 

m 



r(-^)r(i)r(i)r(^(i-.)) 



At next order one finds 



m 



nl-s) r[i + A(i_,)] • 

, r(-5)r(-|) r(i)r(-A(i + ,)) 
n-l-s) r(i-^(i + .)) 



(R37a) 



(R37b) 



(F38) 



(F39) 
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One verifies that the residue of the double pole in s coincide, as expected, with the terms appearing in the expansion 
ofBiE,g). 

(The double-well potential.) This is a symmetric case in which Bi = B2 = B and C = 1 [see equations (2.21), 
(2.22), (2.24a) and (2.58)]. Note that according to (F.17), uq = 1/2 for the double-well potential. We recall the 
relations U{q) = g (1 — g), as well as p = and p{u) = \/l — Au vaUd for the double-weU. At leading order, 
one finds (the duphcation formula of the V function has been used) 

/o(.) = -.-^23-£!ilf:^)im. (F.40) 
°^ ' ^ sr(5/2- 2s) 

For s — > 0, one recovers the contribution l/ig. 

Expanding for s ^ 1 and again using the correspondence (F.7), one obtains the contributions to (3.57) 

2E\n{-E) -2E + 2E \n{-g/2) 2 lnr(i - E) - 2E\n{-2/g). (F.41) 



Similarly, for s ^ 2, 



e^i;^ In I — ) + 17 gE'' (F42) 



70 g'^E^ In ( ) + 236 g^E^. (F43) 



and for s — > 3: 



To obtain the contributions to A{E,g), one has to subtract various contributions of order g^ E^. One of these 
originates from the expansion of the term 

B{E,g)\nB{E,g) ~ gb2{E) \n[E {I + {g/ E)b2{E)] 

~ g'^bl{E)/E r-.s'^g'^E'^/E = 9g'^E^ (F.44) 

[see equations (2.23a) and (2.24a)]. The term B{E, g) \nB{E, g), in turn, comes from the expansion (F.2) of the 
F-function in (F.32). The other terms of order g'^E^ cancel against each other. The result for the term of order 
g^E^, which is 236 — 9 = 227, is in agreement with (2.24b). In order to obtain the results presented in chapter 8.4, 
the contributions of order g^ and 5^ are useful. We have from the contribution of s — > 4: 

1155,3^^1n(f)+l^,3^^ (F.45) 



and from s 



45045 4^,, fEq\ 335183 4^5 

g^E'' In — H g^E^ . (F.46) 
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The complete result for the contour integral of the leading WKB term, through the order of g^, is 

i jdzSo{q,g,E) = l- + i^2E\n(^^-2E 



-g {(iE^ln[ ^ ] + 17E 
+g^ ( 70E^ln( ) +236^;' 



, , a. ,'Eq\ 49843 4 

+g'' 1155£;*ln — ^ + — —E^ 



2 12 



4 /45045 fEg\ 335183 ^. 

5 /II 056 741 fEg\ , 969 969 . 



+9^ ^^^E^ln^-f]+—^E 

(515 954137^,^^ |^^^| ^ 22 309 287^,^^ ^ 

7 /2 151 252 675 o, / Eg\ 469 212 586 743 . 

+^ ( 8 ^ '\-y)+ 448 ^ 

214 886 239 425 q f Eg\ 70 860 581 490 397 



I 32 ^ 1^ r 2688 ^ 1 ' ^^-''^ 



where we neglect terms of order g^ and higher. At next WKB order 
his) can be rewritten as 

72(s) = -5«+i23-i£!tfl£hlli±lf. (F.49) 
^ ^ ^ r(-i -2s) 1 + s ' ^ 



For s — > — 1, one obtains a singular term contributing to the asymptotic expansion of the F-function: 

1 1 
12(s + l) ^ ~12E 



From the expansion for s — > 0, one infers 

■ 2s2 

This yields a contribution to (3.57): 



--^ 1 l-slnl - ) + iis) . (F.51) 



1 , /SflX 11 

:5lnHr +7r.9- (F-52) 



2" V 2 y 6 

To obtain a contribution to A, one must now take into account the correction coming from replacing by -B in the 
expansion (F.32) of r(i - B). At tiiis order only -1/(12B) contributes: 

1 1 9 

+ -,+■■■■ (F.53) 



12B 12i; 4 

It follows tiiat the contribution to A{E, 5) is (1 1/6) - 1/4 = 19/12. Thus, the expansion (2.24b) has been verified. 
For s ^ 1, one infers the contribution to (3.57): 
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Now three terms in (F.32) contribute, involving B up to order 51^ and (equation (2.24a)) one finds 605/12 — 3/2 — 
35/12 + 9/12 = 187/4, in agreement with (2.24b). 

The complete result for the contour integral of the WKB term ^2, up to the order g^, is 



12E ■^U^"l¥)+^ 



+g^ (^Elnf^U^^ 



2 V 2 y 12 

3 l7i^j,2,fEg\ , 4522 ^2 



4 / 45 045 j^3,fEg\ , 743 439 



/2807805 /g,\ 69 706 241 

^ \ 8 \ 2 J 48 

R /88267179 f Eg\ 1097349517 ^ 

7 /2 788 660 875 . f Eg\ 69402 310 265 . 

+^ ( 8 ^ J + 48 ^ 

(^353522522 925 E^ln^^^ 82167 014 713 033 ^ 

where again terms of order g^ and higher are neglected. This result is required for the calculations presented in 
chapter 8.4. 

Without further calculational details, we also give here the contour integral of g^ S4, up to the order g^: 



^^dz[g^Si{q,g,E)] 



Ug ^1015' 



lUOE^ 480^2 480E 



, ,175 , /'Eg\ 17473 



16 V 2 / 288 
+g^ r^Elnf^V^^ 



32 V 2 / 128 

. /1 924 923 o, fEg\ 544644431 , 
+9 ^7; -E^ In Hr + 77:^?^ E 



32 V 2 y 1920 

e / 100 553 453 3 / Eg\ 83 125 560 313 3 



+0° E-'lni — ] + E' 

^ V 32 V 2 y 5760 

7 / 9 526 065 549 ^4 , f Eg\ 645 1 15 861 327 ^3 

+9 ^. E in — H — E 



64 V 2 / 960 



o /1 691 601 686 775 f Eg\ 60 366 482 211 337 .\ 

+g^ E^lnl — ] -\ . (F.56) 

^ V 256 V 2 y 2048 J ^ ' 
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The contour integral of Sq, up to 0{g^), reads: 

1 

9 



jdz[g^S^{q,g,E)] 

__ 31 87 g 359 .g2 15. 2515 

~ ~ 20160 £5 + 4480 £4 + 40320 ^ 128 E"^ ^ 480£ 
_^^5 /159 159 ^ 59 665 081 



,128 V 2 y 7680 

e p5 746 721 25 285 094 891 

^ V 128 V 2 / 23040 

7 / 2 506 538 463 o , / \ 5 392 814 329 807 
+5 TK-o In + E^^HT^ 



128 V 2 / 53760 



o / 758 382 964 625 0, f Eq\ 1884561008165 335 .\ 

E^\d.[ — \ ^ E^ . (¥.51) 

^ V 512 V 2 y 258048 / ^ ' 

The term g^ Ss, upon contour integration, yields to the order g^, 



^j^dz[g'Ss{q,9,E)] = 



127 7381(7 ^ 217 g"^ 



107520 322560 9216 



3377 c/3 ^ 2195/ _ 593 329 5^ ^ 3344883 / 



215040 12288 12288^2 20480 

7 /692049787, f Eg\ 5 827 886 716 943 



2048 V 2 / 2580480 

8 /663 834 081625 ^, f Eg\ 9 724 807 577 177 167 \ 

Eln — H . (F.58) 

^ V 8192 V 2 y 20643840 / 

The leading terms in 1 /E reproduce the coefficients of the asymptotic expansion of the F-function in equation 
(F.2). 

A similar calculation can immediately be repeated for the cosine potential. Note that for the radial Schrodinger 
equation, similar calculations are possible with slight technical modifications, because the WKB expansion is 
vahd for both E and u + 21 large. 



F.7 Alternative Methods for the Contour Integrals 

In chapter F.6, we have considered the evaluation of contour integrals of successive terms in the WKB expansion 
by the method of the Mellin transform. In the current chapter, two alternative methods will be discussed: (i) a 
"subtraction procedure" and (ii) an explicit direct evaluation of the contour integral which is inspired by Lamb- 
shift calculations [31]. 

(Method (i).) We consider the double-well problem and use the convention U{q) = q{l — q). Starting from the 
expansion (B.18), 

2gE Y r(n-l) 



the contour integral of So{q, g, E) can be written as 

- (fdqSo{q,g,E) = - [' dqSo{q,g,E) = - f dq ^U^{q) - 2gE 
9 J 9 Jo 9 Jo 

= - f dq^Jq^{l-qf-2gE 
9 Jo 

^2^ r(n-|) /'dggi-2"(l-g)i-2". (F.60) 
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Here, we have interchanged the infinite sum with the integration, which is not a mathematically rigorous procedure 
in the current situation. Moreover, formally, the resulting integrals do not converge for general positive integer 
n, because of the asymptotic behaviour of the integrand near g ^ and q ^ 1. However, it is easy to evaluate 
the integral in terms of the Beta function for arbitrary n, and to perform an expansion about integer n, by setting 
n ^ n + T]. The following series results, 

^ ^dz So{q, g,E) = ^ + l^2Eln (^^^ - 2E 

+g (^6E'i^\n(^^^+l + Oir,)^ + l7E^ 

+g' (^70E^ jln (^^^ + i + 0(,7)| + 236E^ 



V (ll55i;^{ln(^:^)+i+0 



W + 



49843 . 



+5^ 



/ 45045. 



E9\ , 1 



-E' 
12 

335183 



(F.61) 



When employing the prescription of neglecting the divergent terms in r/, the result (F.47) is recovered. This strategy 
is in formal analogy to the "minimal subtraction" procedure — however, we here do not change the dimensionahty 
of the integration. 

{Method ( ii).) We now discuss the applicabiUty of the so-called e-method for the evaluation of the contour integrals. 
This method has been used for the evaluation of integrals which give rise to logarithmic terms in addition to the 
usual power terms in an asymptotic expansion (for a discussion, including a number of illustrative examples, and 
further references, the reader may consider [31, 32J). We observe that the logarithms originate because of the 
integration regions near q = Q and near g = 1 in (F.60). The divergences are cut off naturally at the scale q~gE, 
but actual divergences result when the expansion in powers of gE is performed. We may therefore consider to 
introduce an arbitrary separation parameter e which separates the regions (0, e) and (1 — e, 1) on the one side from 
the region (e, 1 — e) on the other side. In the latter region, we may safely expand in powers of gE, whereas in 
the first region, we may expand in q (or 1 — q) which is effectively an expansion m gE (this is in analogy to the 
example considered in appendix A of [31]). The arbitrary parameter e cancels if the results of the integrations are 
expanded first in g E, then in e. 

The integration regions (e, 1 — e) gives the following contribution: 



- 'dqVq' {l-qf-2gE 
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(F.62) 
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We should add the result from the regions (0, e) and (1 — e, 1): 



2E I In I - 1 - 2 In e 

+SB" |61n I ^ ) +10- 12 lne + - + ^ 



2 3 I ,/ Jig \ 883 , 70 15 10 1 
+g B'^ 70 In — H 140 In e H h ^ H ^ H 7 

3 4r f E9\ 7934 1155 525 70 35 7 5 

1155 ln^ +_ -2310 lne+ — + — + - + — + — + — 



+g* I 45045 In I — ) + - 45045 Ine + + 



Eg\ 862987 45045 21021 



3003 3465 231 105 9 7 , 

+ ^ + ir-T + 77-^ + T-TT + 77-7 + ^77"^+ ■ (F.63) 



2/16 2£ 4e- 

465 231 105 9 7 
2 e3 ' 8^ 275 ITe 277 16 €» 

Adding the contributions, the dependence on e cancels, and the result (F.47) is recovered. 



F.8 Four-Instanton Coefficients 

We give here the analytic expression for the four-instanton shift of the energy eigenvalues of the double-well oscillator up to eighth order in g. There are four 
infinite series in g, each mixltiplying a specific power of the logarithmic factor A(g) = ln(— 2/g). The shift B^*' (g) is the sum 

fc = 

with the Lk{g) given below. The leading term, for small g, is of the order i*(g)X^(g), 

, , , ,4, x,3, s f 8 352 1156 2 133394 3 16276129 4 550613176 5 

Ls{g) = Sig)X{g)!^--—g+—g-^^g--^g-^^g 

123969635693 „ 40904182753031 , 14705773688939005 „ , „ 1 

g g g + 0(g ) > . (F.65) 

19683 236196 2834352 ^ v» /j 

^-2(3) is given by 



4. ^ 2. ^ f / 352 \ / 1522 1156 

L2(g) = «^s)A'(g) |87+ (^-80- — 7j 9+ (^— + — - 



7 ) 9 

5350 133394 \ 3 / 3332681 16276129 \ 

+ 1-^ ' +1 m2 i^^j ' 

7850396597 123969635693 

+ 
+ 



5832 6561 



/ 158239171 550613176 \ 5 / 117 
I 243 729 7 ® 

/ 4249565076178 40904182753031 \ 7 

V 6561 78732 7 ^ 

/ 1717917660985273 14705773688939005 \ g , „J 

+ 7 9+ 0(9 ) > • (F.66) 

V 78732 944784 r j y yy , 
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The coefficients entering into Li (g) are more complex, 

Li(g) = |*(s)A(s)|(87' + 2C(2)) +(^-34- 1607-^7'■ 
3044 1156 2 289 , A 2 



/ 2564 
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+ 7 7 C(2 9 

V 648 81 486 1944 ^ 

( 2623137449 316478342 550613176 2 137653294 5 

+ 7 7 C(2) g 

\ 9720 243 729 729 •'^ ' ) " 

/ 824874353753 117850396597 123969635693 2 123969635693 \ g 

V 58320 2916 6561 ''' 26244 ^ ') ^ 
/ 2215465381303411 8499130152356 

V 3674160 6561 ''' 

40904182753031 , 40904182753031 , A 7 

7 C(2) g 

78732 314928 -'^ > j " 

/ 1079940312204816343 1717917660985273 

V 44089920 39366 

14705773688939005 2 14705773688939005 , A g „, a 

7 C(2) g + 0(g-^ 
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Tlie series Lo(g) is free of any logaritiims, 
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F.9 Higher-Order Coefficients 

Some instanton coefficients have already been presented in chapters 8.4, 8.7 as well as in appendix F.8. Here, 
we present numerical data for all coefficients up to eight-instanton order, and up to seventh order in g. The ana- 
lytic expressions become rather involved; numerical data in table F.l exhibits the rapid (factorial) growth of the 
coefficients in higher orders in g, as well as the rapid growth for fixed order in g in higher instanton-order. 
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3X10^ 


-1 


4x10'' 


-4 


9X10* 


-1 


8X10^° 


6 


3 


6 


6X10^ 


— 2 


8X10^ 


3.2X10* 


— 1 


0X10^ 


— 2 


9X10^ 


-1 


0x10'' 


-3 


5X10* 


-1 


2X10-*° 


6 


4 


3 


1X10^ 


— 1 


2X10^ 


8.9X10^ 


— 2 


1X10* 


— 1 


6X10^ 


-4 


7x10" 


-1 


4X10* 


-4 


5X10° 


6 


5 


1 


1X10^ 


— 2 


2X10^ 


9.1X10^ 


— 2 


7X10^ 


-3 


8X10* 


-9 


2X10^ 


-2 


4x10'' 


-7.0X10* 


7 





5 


3X10^ 


-3 


7X10^ 


9.1X10* 


-9 


3x10" 


3 


sxio" 


-1 


oxio' 


-3 


0X10* 


-1 


2X10^° 


7 


1 


1 


9X10^ 


— 1 


2x10" 


2.7X10^ 


— 2 


4X10^ 


6 


4X10^ 


-2 


9x10'' 


-1 


0X10^ 


-4 


ixioi° 


7 


2 


3 


1X10^ 


— 1 


9x10" 


3.6X10^ 


— 2 


7X10^ 


4.5X10^ 


-4 


3x10'' 


-1 


6X10^ 


-6 


2X10^° 


7 


3 


3 


0X10^ 


-1 


exio'* 


2.7X10^ 


-1 


6X10^ 


1 


0X10^ 


-4 


1X10^ 


-1 


5X10^ 


-5 


5X10l0 


7 


4 


2 


0X10^ 


-9 


3X10^ 


1.2X10^ 


-5 


3X10^ 


-3 


5X10^ 


-2 


6X10'' 


-8 


8X10* 


-3 


ixioi° 


7 


5 


8 


1X10^ 


-3 


2X10^ 


3.1X10* 


-9 


4X10* 


-3 


1X10^ 


-1 


0x10'' 


-3 


1X10* 


-1 


0X10*° 


7 


6 


2 


3X10^ 


-5 


5X10^ 


3.0X10^ 


-8 


5X10^ 


-7 


4X10* 


-1 


9x10" 


-5 


0x10'' 


-1 


5X10° 


8 





1 


5X10^ 


-1 


2 XIO"* 


3-6X10^ 


-4 


7x10" 


2 


7X10' 


-7 


3X10' 


-6 


6X10" 


-2 


9x10^" 


8 


1 


6 


2X10^ 


-4 


7X10* 


1-3X10^ 


-1 


5x10'' 


7 


1 xio'' 


-1 


8X10^ 


-2 


7X10^ 


-1 


1 XIO** 


8 


2 


1 


2X10^ 


-8 


2X10"^ 


2-0X10^ 


-2 


1 xio'' 


8 


oxio''' 


-2 


1 XIO^ 


-5 


0X10° 


-2 


0X10** 


8 


3 


1 


4X10-^ 


-8 


8X10'* 


1-9X10^ 


-1 


7X10'' 


4 


7X10'' 


-1 


8X10^ 




7X10° 


-2 


2X10** 


8 


4 


1 


IXLO^ 


-6 


2 XIO"* 


1-2 Xio'^ 


-8 


2 Xio'^ 


1 


5x10'' 


-1 


2X10^ 


-4 


3X10° 


-1 


6X10** 


8 


5 


5 


9X10^ 


-2 


9X10* 


4.5X10^ 


-2 


4X10^ 


1 


7X10^ 


-6 


4x10'' 


-2 


2X10° 


-7 


7X10*° 


8 


6 


2 


1X10^ 


-8 


9X10^ 


1.0X10^ 


-3 


9X10^ 


— 4 


0X10^ 


-2 


2x10'' 


-7 


0X10* 


-2 


3X10*° 


8 


7 


5 


2X10^ 


— 1 


4X10^ 


9.7X10^ 


-3 


0X10* 


— 1 


4X10^ 


-3 


9x10" 


-1 


1X10* 


-3 


2X10° 



Table F.l: The e-coefficients (instanton coefficients) determine the resurgent expansion (2.13) 
for the energy eigenvalues of the quantum mechanical double-well oscillator. Here, all coeffi- 
cients up to eighth order in the instanton interaction (up to the "eight-instanton order") and up 
to seventh order in the coupling constant g are considered. Calculations are carried out for the 
ground state. We recall that eo,nki is the coefficient multiplying the term S^{g) X{g) gK The coef- 
ficient eo,n(„_i)o multiphes leading term in the n-instanton order. There is a rapid growth of the 
absolute magnitude of the coefficients in higher orders in g as well as in higher orders in n. 
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